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1. INTRODUCTION 


This is the initial release of a theoretical manual for Version 1.10 of the NASA Rolling Element 
Bearing ANalvsis System (REBANS) family of computer programs. Compared to other bearing 
analysis codes, REBANS offers improved capability to determine the quasi-static response to 
external loads or displacements for three types of high-speed rolling element bearings: 

• angular contact ball bearings (single ball row), 

• duplex angular contact ball bearings (dual ball row), and 

• cylindrical roller bearings, 

by including the effects of bearing ring and support structure flexibility. A finite element 
representation, prepared using Version 4.4 of ANSYS, is used for structure modeling. The 
current analysis system is composed of two main programs, the operation of which are referred to 
separately in a companion user’s manual, are briefly described as follows: 

PREBAN - PREprocessor for Bearing ANalvsis. used to create and modify the two 
necessary input files for FEREBA the main analysis program. This 

interactive code provides extensive on-line help for all model definition 

commands and examines input data for validity. 

FEREBA - Flexibility Enhanced Rolling Element Bearing Analysis, the main analysis 
code, reads the files prepared by PREBAN and performs the requested 
analysis. The program may be run either interactively (with command line 
arguments) or batch, and prints a complete output file of results from the 
analysis. 

The primary focus of this technical manual is the description of the theory that enables the 

flexibility enhanced quasi-static analysis. As such, all discussion in this manual is relative to the 

main analysis code FEREBA a complete operating description of PREBAN is given in the User’s 
Manual. 
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1.1 Development History 

The genesis of FEREBA is taken from one the more recent computer programs developed to 
perform rolling element bearing quasi-static analysis - SHABERTH [1] This code was chosen 
as the development platform for several reasons: 

• modular organization of source, 

• availability of program documentation, 

• use of standard bearing kinematic relations, and 

• ability to analyze the types of bearings required. 

SHABERTH, which is an acronym for SHaft and BEaRmg THermal analysis, has a long history 
of development, starting from the theory originally described by Jones [2], and subsequently 
modified by Harris [3], In addition to the fundamental bearing mechanics analysis, SHABERTH 
also provides the capability for determining the load and temperature distribution for a shaft and 
bearing system, including the effects of lubrication and friction. These additional analysis features 
were not part of the development scope of FEREBA. 

Although SHABERTH provides a sophisticated treatment of a shaft-bearing system, the only 
flexibility in the bearing analysis is due to contact (Hertzian) deformation. For bearings that are 
mounted into rigid housings, without any clearance radially or axially, considering only contact 
deformation is probably reasonable. However, seldom is a bearing installed with such an 
arrangement, and thus the question is raised regarding the influence of clearance and flexible rings 
on the load-deflection characteristics. 

Recognizing the limitations of the rigid-ring theory, several authors have developed extended 
analyses that include the effect of ring flexibility. Filetti and Rumbarger [4] modeled the outer 
ring of a roller bearing with what is essentially curved-beam finite elements. Their calculations 
and experiments showed marked differences in roller load distribution between rigid and flexible 
ring models. Davis and Vallance [5], in a precursor effort to REBANS, developed a method 
using a finite element representation of the outer ring. Their results, for axial preload conditions 
only, showed increased axial motion of the inner ring, compared with rigid-ring calculations. 

REBANS extends these previous efforts to consider general race and housing flexibility, including 
such effects as dead-band and preload springs. The intent of the analysis is to determine the load- 
deflection relationship for system comprised of shafting, the inner and outer rings, the ro ing 
elements, and any significant supporting structure, all represented with 3D finite element models 
In a sense, this system can be thought of as a “control volume for bearing analysis. 


* Numbers in brackets refer to the reference list given in Chapter 9 of this manual 
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1.2 Nomenclature 


Throughout this theoretical description, the following symbols are used to describe physical 
quantities. A brief definition of the terms is also provided. Where possible, common notation 
used in traditional bearing mechanics has been retained. 

A/j Axial distance between ORCC and IRCC under load at bally 

A 2j Radial distance between ORCC and IRCC under load at bally 

B Total curvature for ball bearings 

BD Initial distance between IRCC and ORCC in the unloaded mounted configuration 

CD k Crown drop (including roller crown and raceway crown) of k* slice of all rollers 
measured normal to roller 

CSi Global coordinate system fixed in inertial space with origin at the center of an 

unloaded bearing, with Cartesian axes of X-Y-Z, cylindrical axes of X-R-& 

CS 2 Primary bearing coordinate system which moves with the inner ring and shaft relative 

to the global inertial coordinate system ( CSi ) with five degrees of freedom (Ax, Ay, Az, 
By, 6b), with Cartesian axes of x-y-z, cylindrical axes of x-r-<j> 

D Unloaded rolling element diameter 

D m Diametral distance between rolling element centers in unmounted configuration (pitch 

diameter) 

Eb, E w Elastic (young's) modulus for ball and raceway, respectively 

F C j Centrifugal force on the f ball directed along the positive R axis 

Component of F CJ directed along the x-axis 
Fcrj Component of F CJ directed along the r-axis 

Fiij , Fly Radial force applied to left, right node in j* azimuthal plane of inner raceway in inner 
ring elastic model 

Foij , For/ Radial force applied to left, right node in j* azimuthal plane of outer raceway in outer 
ring elastic model 

F P re Specified total axial preload force 

G Vector of gap separation distances (either axial or radial dimension) 

IRCC Curvature center of a ball bearing inner raceway groove 

Kd Jones axial deflection constant 

Ky , K oj Hertzian contact stiffness between the f 1 ball or roller and the inner, outer raceway 
F , L 0 Effective roller length at inner, outer roller bearing raceway contact 
My , M 0 j Total tilt moment about centroid of j* roller due to local roller contact forces at inner, 
outer raceway 

n Number of rolling elements in a single, primary or duplex bearing 

ORCC Curvature center of a ball bearing outer raceway groove 
p v ( p oJ Total radial force applied to centroid of)* roller by inner, outer raceway 
P d Initial diametral clearance 
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Nomenclature, continued 

Qikj , Qokj Local radial force applied to k* slice of j* roller by inner, outer raceway 

r, , r 0 Radius of curvature of a ball bearing inner, outer raceway grooves 

Rcj Radial component of constraint relation that determines if inner race is unloaded 

Ri , Ro In the initially unloaded condition, locus radius of IRCC s, ORCC s for a ball bearing, 
radial distance from bearing center to inner, outer raceway for a roller bearing 

Xoi , X& Distance along roller axis from roller centroid to center of k** 1 slice measured along 
inner, outer raceway 

X C j Axial component of constraint relation that determines if inner race is unloaded 

X C d Inertial axial distance from the CSj origin to the ORCC of a duplex bearing in the 

unloaded condition (the duplex bearing is initially located a positive distance X cd along 
the X axis relative to the primary bearing) 

Xij Axial distance between ball center and ORCC under load 

X2j Radial distance between ball center and ORCC under load 

a 0 Initial (positive) contact angle for the primary bearing in the unloaded mounted 

configuration defined in CS2 

ay, Ooj Inner, outer race loaded contact angle at the j th azimuthal ball location defined in CS 2 

a p Initial (positive) contact angle for the primary bearing in the unloaded mounted 

configuration, defined in CS 2 , under the action of specified axial preload 

Sdy , Sd 0 j Hertzian contact deflection of the ball and inner and outer race at the j th azimuthal ball 
location defined in CS2 

dn Hertzian contact deflection along the line of contact due to axial preload 

dry Radial elastic displacement of the IRCC at the f ball location defined in CS 2 

dr P i , drpo Radial elastic displacement of inner, outer ring due to axial preload averaged around 
all ball locations 

Sr 0 j Radial elastic displacement of the ORCC at the j 1 * 1 ball location defined in CS/ 

dxy Axial elastic displacement of the IRCC at the j* ball location defined in CS2 

Sx 0 j Axial elastic displacement of the ORCC at the j** 1 ball location defined in CSt 

Sx pi , Sxpo Axial elastic displacement of inner, outer ring due to axial preload averaged around all 
ball locations 

Squy , Sqoig Radial interference between inner, outer contact point at the k* slice on the j* roller 
and radially adjacent point on inner, outer raceway 

duy , Su 0 j Elastic radial displacement of inner, outer raceway in j* azimuthal plane 

dry , dr OJ Radial interference between center inner, outer contact point on j* roller and radially 
adjacent point on inner raceway 

Sruy , dr,rj Elastic radial deflection of left, right node in j* azimuthal plane of inner raceway in 
roller bearing inner ring elastic model 

dr 0 ij , dr OTj Elastic radial deflection of left, right node in j* azimuthal plane of outer raceway in 
roller bearing outer ring elastic model 



Nomenclature, continued 

&e Inertial axial displacement of the CS 2 origin, relative to the CS, origin, required to 

bring a ball into contact with both races during mounting (also the inertial axial 
distance from the ORCC to the IRCC in the unloaded mounted condition) 


A P d 

Arj 

bUrj 

Ax 

by 

AY 
A z 

AZ 
AOj 
Mb, M* 

G,J > 0OJ 

Oy, Gz 

gy, erz 

subscript d 
subscript/) 


Fixed axial displacement of inner ring due to axial preload 

Change in radial clearance due to temperature effects 

Inertial radial displacement of centroid of j 01 roller in j th azimuthal plane 

Total radial displacement of inner raceway in j* h azimuthal plane including elastic radial 

displacement 

Rigid-body displacement of the CS 2 origin in the A’-direction related to shaft/inner ring 
response to applied forces and moments 

Rigid-body displacement of the CS 2 origin in the 7-direction related to shaft/inner ring 
response to applied forces and moments 

Inertial 7-axis displacement of inner ring and shaft 

Rigid-body displacement of the CS 2 origin in the Z-direction related to shaft/inner ring 
response to applied forces and moments 

Inertial Z-axis displacement of inner ring and shaft 

Total tilt rotation of inner raceway at j" 1 azimuthal plane (elastic plus imposed) 

Poisson's ratio for ball and raceway, respectively 

Inertial tilt rotation of centroid of)* roller in j* azimuthal plane 

Azimuthal angle to the j 1 ** ball or roller location in the quasi-static analysis measured 
from the global inertial 7-axis (the planar kinematics used in the quasi-static analysis 
permits O to be considered identical in CS; and CS 2 ) 

Elastic tilt rotation of inner, outer raceway in j* azimuthal plane 

Rigid-body rotation of the CS 2 axes about the 7- or Z-axis related to shaft/inner ring 
response to applied forces and moments 

Rigid-body rotation of the inner ring and shaft about inertial 7- or Z-axis 

Appended to any symbol related to the primary bearing, the subscript d distinguishes 
the variable as being related to the duplex bearing 

Appended to any symbol, the subscript p indicates the variable is related to axial 
preload calculations 
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2. ASSUMPTIONS AND LIMITATIONS 


The theoretical development described in this manual makes assumptions on the bearing system 
behavior and places limits on the extent of applicability as defined below: 

1. The single ball/roller bearing or duplex angular contact ball bearing being studied has a 
nonrotating carrier with fixed supports and a single rotating shaft. The shaft rotates with 
constant angular velocity as required for a quasi-static analysis. 

2. The shaft interfaces with the inner ring by means of a shrink fit with no clearance. Except 
for bearing systems with hollow shafts, elastic displacements of the inner ring and shaft are 
generally considered to contribute negligibly to bearing kinematics. 

3. Inner ring relative structural displacements are obtained from a linear flexibility matrix 
constrained at the interface between the shaft and inner ring in the axial and tangential 
directions. If the shaft finite-element model is sufficiently detailed, potentially significant 
local radial deformations of a hollow shaft are accurately represented. Tangential elastic 
deflections of the inner ring are considered negligibly small. 

4. Force-displacement relations for the outer ring and carrier are described by a single 
stiffness matrix with the associated assumptions such as small deflections and linear elastic 
material properties. The outer ring is connected to the earner by soft springs joining 
coincident nodes across radial and axial gaps. These soft springs must have stiffnesses 
which guarantee that the spurious forces generated across open gaps will be negligible, 
and that the single stiffness matrix describing the outer ring and carrier can be inverted to 
form a flexibility matrix [CO]. Outer ring “rigid-body” displacements result from 
distortions of these soft springs. 

5. Nonlinear effects due to radial and axial contact between the outer ring and the carrier are 
treated with the flexibility matrix using kinematic constraint relations. Due to chamfering, 
any point on the exterior of the outer ring can contact the carrier either in the radial 
direction or in the axial direction, but not in both directions simultaneously. 

6. All forces generated by the cage and all dissipative forces are neglected for the quasi-static 
analysis. The quasi-static analysis, based on the classical approach of Jones [2] and Harris 
[3], determines elastic forces and deflections with balls located in their evenly spaced 
idealized positions. 

7. No transient thermal effects will be considered. Steady-state thermal effects are 
considered in the quasi-static analysis only to the extent of dimensional change. 

Permitted motions of the inner ring consist of five degree-of-freedom (DOF) rigid-body motion 
combined with elastic structural deflections of the inner ring relative to the shaft motions. Rigid- 
body motion of the inner ring is described by small displacements (Ax, Ay, Az) and small rotations 
(0y, 0z) of a reference frame (CS 2 ) fixed at the shaft center and able to move relative to the fixed 
global inertial reference frame (CSi). These five rigid-body freedoms are considered to result 


2-1 


from general shaft bending and axial displacement. Elastic displacements of the inner ring relative 
to (CS 2 ) are determined from a flexibility matrix [CI\. 

2.1 Solution Methodology 

The SHABERTH quasi-static iteration determines rolling element equilibrium and inner ring 
deflections to match applied forces. To incorporate general structural flexibility into the analysis, 
an iterative solution methodology is used that involves dividing the analysis into three sequential 
steps: 

1. With the outer ring fixed in inertial space and the elastic deflection of the inner ring 
specified, ball forces on the inner and outer races are calculated in terms of assumed 
values of the five rigid-body displacements of the inner ring and shaft. This step follows 
essentially the classical quasi-static bearing analysis of Jones [2] and Harris [3] except that 
all displacements, except for the Hertzian ball-race contact deflections, are considered 
fixed during each iteration. The nonlinear equations are solved using a Newton-Raphson 
technique with the required partial derivatives expressed numerically. 

2. Calculated ball forces on the inner race are applied to the flexibility matrix of the inner ring 
and shaft [C7] to determine the corresponding elastic deflections of the inner race 
curvature centers. These ball forces on the inner race are also summed to determine total 
forces and moments, at the origin of the global inertial reference frame (CS 1 ), for 
comparison in Step 3 with the input applied forces and moments. Calculated ball forces 
on the outer race are applied to the flexibility matrix of the outer ring and carrier [C<9], 
using a nonlinear iterative procedure to identify points of contact. The result of this 
iterative procedure is the “rigid-body” and elastic deflections of the outer race curvature 
centers corresponding to the forces applied to the outer race. 

3. Given the deflections of the inner and outer rings corresponding to the ball forces applied 
to the inner and outer races, the differences between the input applied forces and moments 
and the previously calculated forces and moments are used to revise the values of the five 
rigid-body displacements of the inner ring and shaft. This step follows essentially the 
tangent-compliance method of Davis and Vallance [5] except that the system compliance 
matrix does not contain the outer ring flexibilities. Increments of the five displacements of 
the inner ring and shaft corresponding to outer ring displacements are determined from 
“rigid-body” displacements of the outer ring in terms of a least-squares fit. 

Examining the flexibility enhanced iteration, as described by these three steps, at least one and 
possibly two more iterative loops are required to obtain a solution. The outermost loop is 
described by Step 3. If a nonlinear outer race model is used, another loop is required to 
determine displacements with gaps. Thus a total of 4 nested iterations are performed by 
FEREBA. Note that the elastic deflections of the inner ring and outer ring lag the calculated ball- 
race forces by one iteration. If a reasonable error is given to the elastic displacement iteration, 
this difference is minor. 
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2.2 Kinematic Relations 


The kinematics of static contact between a rolling element and race require that, for a ball bearing 
each race curvature center be on a straight line extending from the ball-race contact point through 
the ball center. A similar relationship holds for a roller bearing. For small displacements m he 
tangential (azimuth) direction, a planar relation of both contact forces directed through the ro ling 
element center must hold. This planar kinematic relation is basic to the classical quasi-static 
analysis, and is illustrated in Figure 2-1 for a ball bearing. 


IRCC (loaded) 



Figure 2-1 Planar kinematic diagram for ball bearings 

Therefore, there are no first-order effects on kinematics of small motions in the tangential 
direction There are, however, second-order effects by which tangential displacements of t 
inner ring affect slightly the radial and axial locations of the outer race curvature center in inertial 
space when determined from the inner ring. The radial incremental displacement due to the inner 
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The basic approach of Jones [2] and Harris [3] can therefore be followed with regard to the two- 
dimensional kinematic relations expressed in their planar diagram. However, for the present 
quasi-static analysis, several additional small terms involving rigid-body displacements and 
structural deflections must be included among the primary &c and Sr displacement effects. 

2.3 Use of Finite Element Models 

Flexibility effects of the inner ring and shaft and of the outer ring and earner are represented by 
structural finite-element models (FEMs) generated by Version 4.4 of ANSYS [5], The FEMs are 
input to the bearing mechanics code in the form of stiffness matrices in cylindrical coordinates. 
The stiffness matrices are manipulated as required and inverted in the bearing mechanics code to 
form the desired flexibility matrices in cylindrical coordinates. 

The elastic deflections of the inner ring are determined from a flexibility matrix [C7] in which 
points at the ring/shaft interface are constrained in the tangential and axial directions. The 
resulting calculations for inner ring deflections due to ball-race contact forces are straightforward 
and linear. Assuming a shrink fit without clearance for the ring/shaft interface precludes any 
complicated nonlinear effects for this structural model. 

The elastic deflections of the outer ring are determined from a flexibility matrix [CO] of the outer 
ring and carrier constrained at the bearing supports. Three different models are considered for the 
radial and axial interfaces between the outer ring and carrier: 

A linear model has the outer ring rigidly 
connected to the carrier at the radial and 
axial interface nodes, as shown in Figure 
2-2 for a ball bearing. The outer ring 
undergoes elastic motions influenced by 
the carrier flexibility. 

The first nonlinear model, illustrated in 
Figure 2-3, has the outer ring connected 
to a rigid carrier at the interface nodes by 
soft springs, except for an optional 
preload spring. The outer ring can 
undergo general "rigid-body" and elastic 
motions relative to the rigid carrier. 

The general nonlinear model, as shown in 
Figure 2-4, has the outer ring connected 
to the flexible carrier at the interface 
nodes by soft springs, except for an 
optional preload spring. The outer ring 
then undergoes general "rigid-body" and 
elastic motions influenced by the carrier 
flexibility. 





Ly7\J 



Figure 2-2 Outer Ring/Carrier Model 1 

2 . 



Figure 2-3 Outer Ring/Carrier Model 2 



Figure 2-4 Outer Ring/Carrier Model 3 
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For all three models, the outer ring and carrier FEMs are modeled with separate nodes across the 
gaps at the interfaces. The first model has the outer ring and carrier FEMs rigidly connected (i.e. 
grounded) at the interface nodes. The second model has the outer ring FEM connected to ground 
by soft springs at the interfaces, except for the use of an optional preload spring. The third model 
has the outer ring and carrier FEMs connected by soft springs at the interface nodes, with the 
possible use of finite stiffness preload spring. Clearances between the outer ring and carrier are 
input to the bearing mechanics code for the two nonlinear models. Gaps are positive clearances, 
and interferences are negative clearances. 

For the general nonlinear structural model, the solution technique, as described in section 1 of this 
chapter, initially identifies contact points in the radial or axial directions between the outer ring 
and carrier. Using the contact points, linear outer-ring deflections are calculated from the ball- 
race contact forces using the [CO] flexibility matrix with appropriate constraint relations using an 
iterative procedure. This general method directly accommodates axial preload effects, centrifugal 
force effects, and . nonlinear effects due to clearance between the outer ring and earner 
(deadband). Interference fits are evaluated by using negative clearance distances. 
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3 . REPRESENTATION OF BEARING FLEXIBILITIES 

The structural flexibility characteristics of the bearing inner ring(s), outer ring(s) and carrier are 
developed using Version 4.4 of the ANSYS structural analysis code [5], A total of nine 
configurations are defined for the outer ring(s) and earner for single ball bearing, duplex ball 
bearing set, and cylindrical roller bearing. A total of three configurations are defined for the inner 
ring(s) and shaft of the same ball and roller bearings. 

Stiffness matrices for the ball bearing configuations have a separate node for each outer race 
curvature center ( ORCC) in the outer ring model and for each inner race curvature center ( IRCC ) 
in the inner ring model. These ORCC and IRCC nodes are connected to the ANSYS structural 
model for the appropriate ring by two 3-D spar elements (STIF8) in each azimuthal plane. This 
modeling concept allows a statically-determinate three-hinged arch to approximate the ORCC and 
IRCC planar kinematics due to elastic distortions of the rings. The ORCC and IRCC elastic 
deflections computed in this manner are then essentially independent .of the local Hertzian contact 
deflections at the ball/raceway contact points. With this approach, the total bearing flexibility is 
properly determined by combining the separate effects of Hertzian contact deflections with 
ORCC/IRCC elastic structural deflections. 

Stiffness matrices for the cylindrical roller bearing configuration have two radial degrees of 
freedom (DOFs) for each rolling element on its lines of contact with each raceway. These two 
radial DOFs are the points of application of the equivalent roller Hertzian line contact forces and 
the points where radial structural deflections are determined. At each azimuthal plane, the two 
radial structural elastic deflections are converted to an elastic radial deflection and an elastic tilt 
angle for subsequent use in kinematic constraint relations and roller contact force calculations. 
To separate the effects of Hertzian line contact deflections from raceway structural deflections, 
multiple adjacent nodes along the contact line near each retained raceway DOF are “coupled” in 
ANSYS to move together in the radial direction. 

For the inner ring(s) and shaft, stiffness matrices are generated in ANSYS using cylindrical 
coordinates having the order (R A 0 A Z A ). These stiffness matrices have appropriate fixed 
constraints imposed in ANSYS at the inner ring/shaft interface to permit the shaft/inner ring 
flexibility influence coefficient matrix to be formed. A set of master nodes are specified and the 
substructure option is run to obtain a reduced matrix. In FEREBA the constrained ANSYS and 
substructured stiffness matrix is then transformed into the FEREBA cylindrical coordinates 
system (X F R F Of) prior to eliminating unnecessary DOFs by another matrix condensation. After 
the stiffness matrix is reduced, it is inverted to form the flexibility matrix [C7] which is then 
partitioned for convenient use in subsequent computations. 

Stiffness matrices for the outer ring(s) and carrier are generated in ANSYS in cylindrical 
coordinates with outer ring(s) unconstrained. The earner stiffness matrix, where applicable, is 
supported with appropriate fixed constraints consistent with the design. After being 
substructured, the total ANSYS stiffness matrix is read into FEREBA and then transformed into 
FEREBA cylindrical coordinates. The stiffness matrix for the outer ring(s) is connected to 
ground or to the carrier by soft springs (3 DOF/node) at each candidate contact point except 
preload spring locations and, where applicable, by axial preload springs at the preload location in 
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each azimuthal plane. The constrained stiffness matrix is then reduced to eliminate unnecessary 
DOFs and inverted to form the flexibility matrix [CO\ This matrix is then separated into four 
partitions, in general, for convenient use in subsequent computations. 

The substructure information written by ANSYS is in binary format, and in general, is ordered 
along the wavefront pattern. The REBANS preprocessor PREBAN reads the binary data, and 
sorts the information in ascending order. This substructure matrix is output from PREBAN in 
ASCII format, using the highest precision available for 32-bit computers. The ASCII format was 
selected to enable the stiffness matrix to be transferred between different computing platforms. 


3.1 Bearing/Support Configurations 

The nine bearing/support configurations for the outer ring and earner are designated by IBS COR 
(index) values one through nine. The three bearing/support configurations for the inner ring and 
shaft are designated by IBSCIR values one through three. In the following figures which define 
these 12 configurations, master DOFs in both ANSYS and FEREBA are listed with n referring to 
the number of rolling element azimuthal planes. Node ordering and, where applicable, gap 
ordering in each azimuthal plane are also indicated. These specified node and gap orderings 
associated with the IBSCOR/IBSCIR indices are used in several subsequent numerical 
computations. 

1. Single ball bearing with flexible outer ring and/or carrier and no deadband (EBSCOR = 1). 
This model is linear. There are no gaps between the outer ring rigid housing, and at each 
azimuthal position, 2 master nodes are required. After reduction in FEREBA, four DOFs 
remain at each ball location (2 axial, 2 radial). 

Active DOFs 

ANSYS 6-n (X, Y, Z) 

FEREBA 4 n (X, R) 

Active Nodes fDOFs) 

X: 1,2 (1,4) 

R: 1,2 (2, 5) 

Reduced DOFs 
3*, 6 

* DOF 3 is deleted 
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///////, 


2. Single ball bearing, flexible outer ring with deadband, carrier assumed rigid (IBSCOR = 2). 
This model is nonlinear due to deadband contact, which is dependent on ball loads applied to 
the outer race. There are four gaps between the outer ring and rigid housing, and at each 
azimuthal position, 6 master nodes are required. After reduction in FEREBA, eight DOFs 
remain at each ball location (4 axial, 4 radial). 



Active DOFs 

ANSYS 18-«(X, Y,Z) 
FEREBA 8 •« (X, R) 

Active Nodes (DOFs) 

X: 1,2,3,6(1,4,7,16) 

R: 1,2,4,5(2,5,11,14) 

Reduced DOFs 

3*, 6, 8, 9, 10, 12, 13, 15, 17, 18 
* DOF 3 is deleted 


3. Single ball bearing with flexible outer ring and carrier with deadband (IBSCOR = 3). This 
model is nonlinear due to deadband contact, which is dependent on ball loads applied to the 
outer race, and also due to the preload spring (nodes 6 - 10), which can bottom. There are 
four gaps between the outer ring and carrier, and at each azimuthal position, 10 master nodes 
are required (6 on the outer ring, 4 on the carrier). After reduction in FEREBA, twelve DOFs 
remain at each ball location (6 axial, 6 radial). 



Active DOFs 

ANSYS 30 n (X, Y, Z) 

FEREBA 12'W (X, R) 

Active Nodes (DOFs ) 

X: 1,2,3,6,7,10 
(1,4, 7, 16, 19, 28) 

R: 1,2, 4, 5, 8, 9 

(2, 5,11, 14, 23, 26) 

Reduced DOFs 

3*, 6, 8, 9, 10,12, 13,15, 
17, 18, 20,21,22, 24, 25, 
27, 29, 30 
* DOF 3 is deleted 
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4. Duplex ball bearing set with flexible outer rings and carrier with no deadband (IBSCOR = 4). 
The analysis is restricted to having n equal in both bearings. This model is nonlinear because 
the preload spring (3 - 6) could bottom. There is only one gap between the outer rings, with 
the primary bearing free to move axially (but not radially). At each azimuthal position, 6 
master nodes are required (3 for each ring). After reduction in FEREBA, ten DOFs remain at 
each ball location (6 axial, 4 radial). 


Active DOFs 

ANSYS 18-h (X, Y, Z) 

FEREBA 10-w (X, R) 

Active Nodes fDOFs) 

X: 1-6(1,4,7,10,13,16) 

R: 1,2,4,5(2,5,11,14) 

Reduced DOFs 
3*, 6, 8, 9, 12*, 15, 17, 18 
* DOFs 3 & 12 are deleted 
Note the relative positions of “primary” and “ duplex ” bearings. 

5. Duplex ball bearing set with flexible outer rings with deadband and with carrier assumed rigid 
(IBSCOR = 5). The analysis is restricted to having n equal in both bearings. This model is 
nonlinear due to deadband contact and the possibility that the preload spring (6 - 9) could 
bottom. There are seven gaps between the outer ring and rigid housing, and at each azimuthal 
position, 12 master nodes are required (6 on each bearing). After reduction in FEREBA, 
sixteen DOFs remain at each ball location (8 axial, 8 radial). 



Active DOFs 



The relative positions of the “ primary ” and “ duplex ” bearings is 


ANSYS 3 6-n (X, Y, Z) 

FEREBA 16 « (X, R) 

Active Nodes fDOFs l 

X: 1, 2, 3, 6, 7, 8, 9, 12 

(1,4, 7, 16, 19, 22, 25, 34) 

R: 1,2,4,5,7,8,10,11 

(2, 5, 11, 14, 20, 23,29, 32) 

Reduced DOFs 

3*. 6, 8, 9, 10, 12, 13, 15, 17, 
18, 21*. 24, 26, 27, 28, 30, 31, 
33, 35, 36 

* DOFs 3 & 21 are deleted 
the same as configuration 4. 
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6. Duplex ball bearing set with flexible outer rings and earner with deadband (IBSCOR 6). 
The analysis is restricted to having n equal in both bearings. This model is nonlinear due to 
deadband contact and the possibility that the preload spring (6 - 9) could bottom. There are 
seven gaps between the outer ring and earner, and at each azimuthal position, 18 master 
nodes are required (6 on each bearing and 6 on the earner). After reduction in FEREBA, 
twenty two DOFs remain at each ball location (10 axial, 12 radial). 

Active DOFs 

ANSYS 54 •« (X, Y, Z) 

FEREBA 22 n (X, R) 



Active Nodes (DOFs) 

X: 1,2,3,6,7,8,9,12,13,18 
(1, 4, 7, 16, 19, 22, 25, 34, 
37, 52) 

R: 1,2,4,5,7,8,10,11,14, 

15, 16, 17 (2, 5, 11, 14, 
20, 23, 29, 32, 41, 44, 47, 
50) 

Reduced DOFs 

3*, 6, 8, 9, 10, 12, 13, 15, 17, 
18, 21*, 24, 26, 27, 28, 30, 31, 
33, 35, 36, 38, 39, 40, 42, 43, 
45, 46, 48, 49,51,53, 54 


* DOFs 3 & 21 are deleted 

The relative positions of the “primary" and “ duplex ” bearings is the same as configuration 4. 

7. Cylindrical roller bearing with flexible outer ring and earner, no deadband (IBSCOR = 7). 
This model is linear. There are no gaps between the outer ring rigid housing, and at each 
azimuthal position, 2 master nodes are required. After reduction in FEREBA, two radial 
DOFs remain at each roller location. 


Active DOFs 



ANSYS 6-n (X, Y, Z) 
FEREBA 2 n (R) 

Active Nodes fDOFs) 

X: none 
R: 1,2 (2, 5) 

Reduced DOFs 


1 , 3 , 4, 6 

The radial force and corresponding moment about the axial centerline will be represented by 
two unequal radial forces at master nodes 1 & 2. 
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8. Cylindrical roller bearing with flexible outer ring with deadband and with carrier assumed rigid 
(ffiSCOR = 8). This model is nonlinear due to deadband contact. There are two gaps 
between the outer ring and rigid housing, and at each azimuthal position, 4 master nodes are 
required. After reduction in FEREBA, four radial DOFs remain at each roller location. 



1 2 


Roller 



Active DOFs 

ANSYS 12-«(X, Y,Z) 
FEREBA 4 n (R) 

Active Nodes (DOFs) 

X: none 

R: 1,2,3,4(2,5,8,11) 

Reduced DOFs 
1,3,4, 6, 7, 9, 10, 12 


9. Cylindrical roller bearing with flexible outer ring and carrier with deadband (IBSCOR = 9). 
This model is nonlinear due to deadband contact. There are two gaps between the outer ring 
and carrier, and at each azimuthal position, 6 master nodes are required (4 on the bearing and 
2 on the carrier). After reduction in FEREBA, six radial DOFs remain at each roller location. 



Active DOFs 

ANSYS 1 8*n (X, Y, Z) 
FEREBA 6 n (R) 

Active Nodes (DOFs) 

X: none 

R: 1-6(2,5,8,11,14,17) 

Reduced DOFs 
1,3,4, 6, 7, 9, 10,12, 

13, 15, 16, 18 
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10. Single ball bearing with flexible or rigid shaft (ffiSCIR = 1). This model is linear. The inner 
ring is assumed to be fixed to the shaft. After reduction in FEREBA, four DOFs remain at 
each ball location (2 axial, 2 radial). 



Active DOFs 

ANSYS 6 77 (X, Y, Z) 

FEREBA 4 7? (X, R) 

Active Nodes (DOFs) 

X: 1,2 (1,4) 

R: 1,2 (2, 5) 

Reduced DOFs 
3*, 6 

* DOF 3 is deleted 


11. Duplex ball bearing set with 
flexible or rigid shaft (IBSCIR = 2). The analysis is restricted to having n equal in both 
bearings. This model is linear, with the inner rings fixed to the shaft. After reduction in 
FEREBA, eight DOFs remain at each ball location (4 axial, 4 radial). 



Active DOFs 

ANSYS 12-w (X, Y, Z) 
FEREBA 8-/7 (X, R) 


Active Nodes fDOFsl 
X: 1, 2, 3, 4 (1, 4, 7, 10) 
R: 1,2,3,4(2,5,8,11) 


Reduced DOFs 
3*, 6, 9*, 12 

* DOFs 3 & 9 are deleted 

Note the relative positions of “primary" and “ duplex ” bearings. 
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12. Cylindrical roller bearing with flexible or rigid shaft (IBSCIR - 3). This model is linear. The 
inner ring is assumed to be fixed to the shaft. 



Active DOFs 

ANSYS 6 n (X, Y, Z) 

FEREBA 2 n (R) 

Active Nodes (DOFs) 

X: none 
R: 1,2 (2, 5) 

Reduced DOFs 
1, 3, 4, 6 


The radial force and corresponding moment about the axial centerline will be represented by 
two unequal radial forces. 


3.2 Inner Ring/Shaft Finite-Element Models 

The ANSYS structural model of the inner ring(s) and shaft is generated in cylindrical coordinates 
having the ANSYS order (R A 0 A Z A ). The shaft/inner ring free-free stiffness matrix must be 
constrained such that the relative elastic deflections from the flexibility matrix [C7] can subse- 
quently be added properly to the rigid-body inertial displacements defined in the shaft coordinate 
system CS 2 . For single ball or roller bearings or primary bearings of a duplex set, the axial and 
tangential DOFs are fixed against displacement at the interface between the shaft and inner rings. 
To represent a solid shaft without a finite-element model of the shaft, the radial DOFs could also 
be constrained. 

If the local radial flexibility effects of a hollow shaft are significant, the shaft must be modeled 
with ANSYS general shell elements (STIF63) or equivalent. If elastic straight pipe elements 
(such as ANSYS STIF16) or equivalent beam elements are used, local radial deformations are not 
available. The shaft would be modeled most accurately with the same radial definition of elements 
as used for the inner ring. This would probably require two to four shell elements between each 
rolling element for each axial section. In the axial direction, the mesh refinement must be 
compatible with the radial definition of elements. 

The master or retained DOFs from the ANSYS inner ring/shaft stiffness matrix must be those 
defined for the specific bearing/support configurations. Unnecessary DOFs not constrained at the 
shaft/inner ring interface or those coupled along the roller/raceway contact line are eliminated in 
ANSYS by static condensation. 

For ball bearings, the inner raceway groove is a circular section for containing the spherical balls. 
The candidate ball contact points within this circular race are defined by nodes corresponding to 
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an assumed constant contact angle of, for example, ±20°. Although the range of typical ball 
contact angles is 0° to ±40°, applying the ball contact force components to a single node is 
sufficiently accurate for determining elastic deflections of the inner rings. By definition, the 
positive contact angle for ball bearings corresponds to a single bearing or to the “primary” bearing 
in a duplex set; the negative contact angle corresponds to the “duplex” bearing in a duplex set. 

With regard to the ball bearing kinematic relations, the most important displacements generated 
by the shaft/inner ring flexibility matrix are those describing the locations of the inner-race 
curvature centers ( IRCCs ). It is therefore necessary to include a node for the IRCC at each ball 
location. The three DOFs for each IRCC node are then connected rigidly to inner race nodes 
adjacent to the ball contact node in the azimuthal plane using ANSYS 3-D spar elements (STIF8). 
Note that the tangential DOF at each IRCC node has no stiffness normal to the azimuthal plane. 

For ball bearings, forces applied to the inner race flexibility matrix [Cl] are calculated in CS 2 
components at each rolling element contact point. The desired elastic deflections are then the 
product of matrix [Cl] times the vector of forces applied to the inner raceways. Matrix [Cl] is 
used only to determine the elastic deflections of the inner ring(s) relative to the CS 2 shaft axes. 
The static force balance for comparison with the applied forces and moments is determined from 
the actual force direction and location and not from the nodal geometry of the flexibility matrix. 

For roller bearings, all forces are defined directly in the inertial radial direction at each azimuthal 
plane. The static force balance for comparison with the applied forces and moments is determined 
using inertial radial directions only. 

In a ball bearing, Hertzian contact deflections at the ball-race contact points are not represented 
by the relatively coarse FEM mesh used to predict elastic deflections of the IRCCs. Using 3-D 
spar elements connecting IRCC nodes with nodes adjacent to the ball-race contact points 
generates structural deflections of the race independently of Hertzian contact deflections. 
Therefore, the Hertzian contact deflections for ball bearings are calculated separately by 
traditional bearing mechanics relations and superimposed on the IRCC elastic displacements for 
the bearing kinematic equations. 

The rolling element contact points for applied forces with cylindrical roller bearings are identical 
to the locations where elastic deflections are required. The radial line contact force and 
misalignment moment between each roller and the inner raceway is approximated by two unequal 
radial forces. These forces are applied at two nodes on the raceway, in each azimuthal plane, 
which are symmetrically located on either side of the raceway centerline a distance equal to 
approximately 30-percent of the roller effective length. This spacing, which must be considered 
in defining the ANSYS mesh for the inner ring model, corresponds to a uniformly loaded beam on 
two supports symmetrically located such that the negative moment at the support equals the 
maximum positive moment between the supports. This spacing is one approach for representing a 
continuous line contact force with two concentrated forces. 

In the roller bearing model, Hertzian line contact deflections are separated from elastic ring 
deflections by means of a special feature in the ANSYS model. Along each line of contact 
between roller and raceway, several nodes adjacent to the retained nodes are constrained radially 
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by the ANSYS “couple” command. By constraining several adjacent nodes to move radially 
together, the Hertzian line contact force is essentially spread along the effective length of the 
roller instead of being concentrated at the two retained nodes. The use of this modeling technique 
is described and illustrated in the User’s Manual. 

3.3 Outer Ring/Carrier Finite-Element Models 

The ANSYS structural model of the outer ring and carrier is generated in cylindrical coordinates 
having the ANSYS order (R A 0 A Z A ). The outer ring(s) are unconstrained, but the carrier support 
DOFs are constrained to zero deflection. The master or retained DOFs from the ANSYS outer 
ring/carrier stiffness matrix must be those defined for the specific bearing/support configurations. 
Unnecessary DOFs not constrained at the earner support DOFs or those coupled along the 
roller/raceway contact line are eliminated in ANSYS by static condensation. 

•For ball bearings, radial and axial gaps are defined in the ANSYS model by coincident but 
separate nodes on the outer ring(s) and carrier. Specified initial gap values for each potential 
contact are subsequently input to FEREBA. Because the outside comers of the outer ring(s) are 
chamfered, the comer points on the outer ring cannot contact the carrier. Radial contact can 
therefore occur only on a radial face, and axial contact can occur only on an axial face. 
Cylindrical roller bearings are defined to have only radial contact. 

In FEREBA, gap closure between the outer ring(s) and earner causes contact forces between the 
two bodies. These contact forces are calculated from radial/axial constraint equations in which 
the radial/axial displacements on both sides of a closed gap are constrained to move together. 
The only nonlinear aspects of this approach are identifying which gaps are closed and nulling any 
impossible tensile forces at the local axial preload springs for ball bearings. Once a set of gap 
closures is assumed, the elastic displacement analysis is linear except for liftoff of the local axial 
preload springs. 

Contact forces between the outer ring and carrier are not calculated in terms of nonlinear Hertzian 
contact force/displacement relations. Such calculations require assumptions of the ideal nature of 
the contact which would not be appropriate for these conditions. In the radial direction, contact 
between convex and concave cylindrical surfaces with essentially equal diameters constitutes 
neither ideal point nor ideal line contact. In the axial direction, contact between two flat surfaces 
ran not be considered Hertzian contact. Including nonlinear Hertzian contact force/displacement 
relations would therefore add unwarranted complexity to the analysis (including potential 
Newton-Raphson convergence problems) without improving the fidelity of the bearing model. 

For ball bearings, the basic purpose of the ANSYS model of the outer ring and carrier is to 
produce elastic axial and radial deflections of each outer race curvature center ( ORCC ). Each 
ORCC point is represented in the ANSYS structural model as a node located properly relative to 
the raceway groove at each azimuthal position. The ORCC node is connected to the structural 
model by two 3-D spar elements (STIF8) in the azimuthal plane as with the ANSYS model of the 
inner ring and shaft. The spar elements connect ORCC nodes with nodes adjacent to the ball-race 
contact points. The resulting elastic deflections are essentially independent of the Hertzian 
contact deflections of the ball-race contact points. 
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As with the inner ring/shaft model for cylindrical roller bearings, the two retained nodes defining 
the contact points between roller and outer raceway are symmetrically located on either side ot 
the centerline a distance equal to about 30-percent of the roller effective length. The nodes 
adjacent to these two retained nodes are coupled in ANSYS so that the radial DOFs move 
together with the retained nodes. Refer to the User’s Manual for more details and examples of 

this modeling technique. 


3.4 Coordinate Transformations 

The stiffness matrices generated in ANSYS must be transformed from the ANSYS cylindrical 
coordinate system {q A } into FEREBA cylindrical coordinates {q F }- The definitions of cylindrical 
and Cartesian coordinates for ANSYS and FEREBA are given in equations (3-1). 


Ra - Rf 
0a = Of 
Z a = X F 


(3-1) 


where R 0, and Z or X are the radial, azimuthal (tangential), and axial coordinates. The 
corresponding DOF orderings are {q A } T - {Ra Oa Z a } and {q F } - {X f Rf ^' nc ® t^eDOF 
order is not the same, a transformation matrix to convert a stiffness matrix from ANSYS to 
FEREBA cylindrical coordinates in matrix notation is given by 


Ra 


0 1 0 

Xf 

Qa 

> = 

0 0 1 

Rf ' 

Xa \ 


.1 0 0. 

.9^ 


(3-2) 


which in terms of the generalized coordinates may be expressed as 

{qA} ~ [T] {qp} 


(3-3) 


The transformation of the ANSYS stiffness matrix into FEREBA coordinates is accomplished by 
the matrix triple product [T] t *[K)*[T\ where matrix [7] is the transformation matrix of equation 
(3-3). This matrix operation is performed in FEREBA by subroutine TRANAF, using an eflficien 
computational procedure to minimize storage requirements. 


3.5 Constraining Free-Free Outer Ring 




Subroutine ADDSTF operates on the unreduced stiffness matrix of the outer ring(s) and carrier in 
FEREBA cylindrical coordinates. The stiffness matrix for the outer ring(s) is connected to 
ground or to the carrier by soft springs (3 DOF/node) at each candidate contact point except 
preload spring locations and, where applicable, by axial preload springs at the preload location in 
each azimuthal plane. The spring stiffnesses are added to the unreduced stiffness matrix in 
locations specified for each of the nine bearing/support configurations. 

The stiffness magnitude of the soft springs is calculated as a factor multiplied by a typical outer 
ring stiffness term. The factor is currently selected as 1.0 x 10" 6 consistent with a similar ANSYS 
approximation. The typical outer ring stiffness term is the radial DOF at the first ball contact 
point for ball bearings and the radial DOF at the second roller contact point for cylindrical roller 
bearings. 

For ball bearings only, the total axial preload spring stiffness input to FEREBA is divided by the 
number of rolling elements to give the magnitude of the individual preload spring stiffness at each 
azimuthal plane. The individual preload springs are uncoupled from one azimuthal plane to 
another. Preload springs are not allowed with a roller bearing. 


3.6 Condensing and Inverting Matrices 

After the stiffness matrices have been transformed to the FEREBA coordinate system, DOF that 
are not used in the subsequent kinematic relationships are removed either by static condensation 
or deletion. The unrequired DOF are listed with the 12 bearing/support combinations described in 
section 3.1. An efficient static condensation procedure has been incorporated into FEREBA 
which reduces DOF one at a time. In partitioned form, the general static condensation problem 
for the stiffness matrix may be stated as 


"[*>] . 1(^)1 
.{AT*} K„ J1 , J 1 F t I 


where x s is the coordinate to be condensed and K ss is the term to be eliminated from the stiffness 
matrix. Note that submatrices {AT,*} and {K JS } are vectors. For convenience, the DOF to be 
condensed is the last term in equation (3-4), practical implementation of this algorithm proceeds 
from the largest DOF to the smallest. 

The second row of (3-4) may be written as 

*, " 4 ~{ f , - {*„>{*,)) C 3 - 5 ) 

^SS 

which expresses the constraint relation between x s and {x y }. Normally, as is the case in FEREBA, 
F, is zero. 
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Expanding the first row of (3-4), and substituting the relation given by (3-5) results in the 
condensed expression for the statics problem expressed in (3-4) as 


.[*,*]- 


{K JS ){K sk } 


I {*;} = {Fj} 


{KjsJF, 

K„ 


(3-6) 


which by inspection, the reduced order stiffness matrix is obtained from (3-6) as 


[k;] = [K jk ] - 


{KjsHK*} 

K ss 


(3-7) 


where [ ] is the condensed matrix. In FEREBA, subroutine MATRED performs this reduction 

by sequentially removing each required DOF and shifting the remaining rows and columns of the 
matrix up and over by one term. 

For some of the ball bearing ring configurations, DOFs are deleted, rather than condensed. 
Mathematically, this procedure applies a fixed constraint to the system, thus only unrequired 
tangential DOF are treated in this manner. In FEREBA, deleted DOFs bypass the reduction 
procedure and are eliminated by row and column shifts. 


Once the stiffness matrices have been reduced, compliance matrices are obtained by using 
standard mathematical subroutines for inversion. FEREBA uses the routine DLFIRG from the 
commercially available IMSL* subroutine library, however, any routine capable of inverting 
symmetric double precision matrices can be used instead. 


3.7 Partitioning Compliance Matrices 

After the stiffness matrices for the inner ring/shaft and outer ring/carrier have been reduced and 
inverted, the resulting flexibility influence coefficient matrices are partitioned for convenience in 
subsequent computations. Subroutine CEP ART extracts influence coefficients from matrix [C7] 
which define IRCC or roller deflections in terms of rolling element forces applied to the inner 
raceway. Similarly subroutine COPART extracts influence coefficients from matrix [CO] required 
for subsequent calculations. 

The four required partitions of matrix [CO] are defined as follows: 

COP 12 = flexibility coefficients defining ORCC or roller deflections in terms of rolling 
element forces applied to the outer raceway; 

COP 13 = flexibility coefficients defining ORCC or roller deflections in terms of 
constraint forces between outer ring and carrier at the candidate contact points; 

COP32 = flexibility coefficients defining deflections at the candidate contact points 
(gaps) in terms of rolling element forces applied to the outer raceway; 

COP33 = flexibility coefficients defining deflections at the candidate contact points 
(gaps) in terms of constraint forces between outer ring and carrier at the 
candidate contact points. 


The IMSL MATH/LIBRARY is available from Visual Numerics, Inc. 
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Subroutines CODIM and CIDIM perform the functions of assigning dimensions for each of the 
one-dimensional and two-dimensional arrays involved with the bearing flexibility computations. 
The array dimensions are determined as the product of the number of rolling elements times 
integers defined for each of the twelve bearing/support configurations. The current version of 
REBANS limits the number of rolling elements to 20 for a single bearing row (40 total for a 
duplex bearing system). 

3.8 Defining Gap Clearance Vector 

For seven of the nine outer ring configurations, a combination of axial and radial gaps are user 
specified by entering values for clearance and diameters. This data is converted to the required 
gap clearances (or interferences in special cases) by subroutine GAPDEF. Currently, this routine 
applies constant radial and/or axial clearance values at all azimuthal positions. More general 
versions can use nonuniform clearance values if desired. 

Depending on the configuration, from two to seven gap clearances require definition. With the 
uniform gap limitation, only four gap values are possible, three axial and one radial. The axial 
gaps correspond to the distance between the outer ring and carrier/housing in the load direction 
for a ball or duplex bearing, denoted as C A , the corresponding clearance on the unloaded side of 
the bearing, set equal to unity (1.0), and the distance between outer rings in a duplex set, referred 
to as Cd ■ The single radial gap is the radial clearance between the outer ring and carrier/housing, 
represented as C R . For the seven IBSCOR configurations which have gaps as numbered in section 
3 . 1 , the clearances are set to the following values: 

IBSCQR = 2 or 3 (ball bearing with deadband) 


I Gap number and direction I 

n— 



DSSIil 


c R 

c R 

c A 1 


IBSCOR = 4 /duplex bearing without deadband) 

Only one gap is required, between the two outer rings, equal to C D 
IBSCOR = 5 or 6 /duplex bearing with deadband) 


1 Gap number and direction | 

1 (axial) 

2 (radial) 

3 (radial) 

gypsEirt! 

5 (radial) 

6 (radial) 

7 (axial) 

1.0 

C R \ 

C R 

C D 


Cr 

C A 


IBSCOR = 8 or 9 (roller bearing with deadband) 


| Gap number and direction f 


BBn 


Cr _J 
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4. COORDINATE SYSTEMS 

This section describes the coordinate systems used to develop the kinematic relations for ball and 
roller bearing analysis. For both bearing types, a fixed inertia system will be defined, from which 
any local systems will be referenced. 

4.1 Angular Contact Ball Bearing Coordinate Systems 

Two separate coordinate systems are defined for the major elements of the ball bearing. One is 
for the inner ring; the other is for the carrier and outer ring. The coordinate system describing the 
outer ring and carrier is the global reference frame fixed in inertial space. The reference frame for 
the inner ring can move relative to the global inertial frame. 

Each coordinate system consists of both right-handed Cartesian (X-Y-Z) axes and cylindrical (X- 
R-0) axes. The origins of each pair of Cartesian and cylindrical axes are coincident, and the X 
axes of each pair are collinear. For all axis systems, positive X is defined in the general direction 
of the applied axial thrust load. The transformations between each Cartesian and corresponding 
cylindrical coordinates are defined as follows: 

Y = R- cos(O) 

Z = i?sin(<D) 

R = Ycos(0) + Z-sin(0) 

<D = tan \Z/Y) 

Note that Xis identical for both Cartesian and cylindrical axes and that 0 is measured from Y with 
the positive sense defined by the rotation of Y into Z. It should be noted that the azimuth angles 
(0) defining ball locations are identical in both the global inertial reference frame (CSi) and the 
shaft/inner ring moving reference frame (CS 2 ). This is because the azimuth angles for the two 
coordinate systems are coincident in defining the initial or reference ball locations before loading. 
Subsequent elastic distortions resulting from loading are measured, for the small displacements 
assumed in bearing mechanics theory, relative to the reference ball locations and do not involve 
any change in the reference azimuth angles. 

4.1.1 Global inertial reference frame (CSi) 

The global reference frame CSi is fixed in inertial space at the center of the unmounted bearing 
and at the center of the carrier bore. X is directed along the unloaded shaft with positive X in the 
direction of the applied axial thrust load. The locus of curvature centers for both inner and outer 
rings is initially located in the Y-Z plane before mounting and loading. 

The initial diametral clearance P d is defined in this reference frame, where P d is total clearance 
across a diameter built into the bearing by the manufacturer. Half of P A is located on either side of 
the origin by definition. 
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The carrier and outer ring are located directly in this inertial reference frame. The carrier support 
points are fixed in inertial space. Any initial clearance between the outer ring and the carrier in 
the unloaded configuration is defined directly in this coordinate system. The general small 
displacements of the outer ring and carrier are determined from a flexibility matrix defined in this 
global inertial reference frame. The locus of outer-race curvature centers is considered part of the 
outer ring. All applied forces and moments, including the centrifugal force on each ball, are 
defined in the global inertial reference frame. 

4.1.2 Shaft/inner ring reference frame ( CS 2 ) 

The origin of CS 2 'is fixed at the center of the loaded shaft during operation. The axial coordinate 
x is directed along the loaded shaft and is positive in the general direction of the applied axial 
thrust load. The origin is initially translated along the X-axis in the positive X direction a distance 
(Ae) equal to half the "free play" of the bearing. This initial position, which has the balls in 
contact with both races, is the starting point for the Jones and Harris bearing analysis. The initial 
mounted contact angles at the inner-race and outer-race contact points are equal and given as 
follows for the unloaded configuration: 

cos(a° ) = (2BD - Pd)/2BD 
Ae = BD sin(a°) 

BD = r 0 + r t - D 

Here r 0 and r, are the radii of curvature for the outer and inner races, respectively, and D is the 
unloaded ball diameter. 

The origin of this shifted moving coordinate system is located in the global inertial system ( X-Y-Z ) 
by three small displacements (Ax+ Ae Ay, A 2 ) and two small angles (0y, 9z). These five DOFs 
(excluding Ae) are the unknown independent variables in the nonlinear equations evaluated by the 
iterative tangent-compliance solution technique. The relationship between these two coordinate 
systems is illustrated in Figure 4-1 . 

The orthogonal transformation matrix from CSi to CS 2 is as follows with the small angle 
assumption: 


X 


1 

A 

e/ 


X 

y 

► — 

-e z 

1 

0 


Y ’ 



A 

0 

1 


z 


This coordinate system locates all points on the inner ring including the locus of points defining 
the inner-race curvature centers which lie initially in the y-z plane. Note that the inner-race 
curvature centers move out of the .y-r plane under general loading due to elastic displacements of 
the shaft/inner ring structure. The Hertzian contact deflections and the ball/race contact angles 
are also defined in CS 2 . This coordinate system is the basic reference frame used by Jones and 
Harris for defining their internal forces and kinematic relations. This coordinate system does not 
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rotate with orbital speed, as it is not fixed to the shaft/inner ring structure. This reference frame 
moves only with the five DOFs defined previously. 



4.1.3 Use of coordinate systems 

Kinematic constraints leading to ball forces applied to inner and outer races are rigorously 
developed from CSi and CS 2 . Ball forces applied to the inner race and the corresponding elastic 
deflections are defined in CS 2 , however, total forces and moments are assumed to be in CS i. This 
simplifies the structural analysis by allowing the torque from ball forces to be zero. This 
approximation thereby allows the bearing system to be defined by the standard five degrees of 
freedom instead of six. 

Ball forces applied to outer race for determining elastic deflections are assumed to be in CSy This 
simplifies the structural analysis by allowing tangential forces to be zero. This approximation 
thereby allows the nonlinear structural analysis to be defined by radial and axial freedoms only. 

These approximations are conventionally used in small-deflection structural analysis. The 
increased kinematic precision used to determine ball forces in CS 2 is considered necessary to 
evaluate the effects of structural flexibility. 
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4.2 Cylindrical Roller Bearing Coordinate Systems 

A single inertial reference frame is used to define radial displacements and tilt rotations in the X-R 
azimuthal planes for all major roller bearing elements. All radial translations and tilt rotations of 
the inner ring and shaft, the outer ring and carrier, and the rollers are defined directly in this single 
inertial reference frame. The inertial reference frame consists of both right-handed Cartesian ( X - 
Y-Z) axes and right-handed cylindrical (X-/?-<J>) axes. The origins of the Cartesian and cylindrical 
axes are coincident, located at the center of the unloaded bearing, with half the diametral 
clearance on each side of the origin. The X axes of both coordinate systems are collinear. Since 
cylindrical roller bearings do not accept axial thrust loads, the positive X direction is defined 
arbitrarily to complete the right-handed coordinate system. The transformations between the 
Cartesian and cylindrical coordinate system are defined as follows: 

Y = R cos(0) 

Z = R sin(<Z>) 

R = Fcos(<J5) + Z sin(<P) 

0 = tan ~\Z/Y) 

Note that X is identical for both Cartesian and cylindrical axes and that O is measured from Y with 
the positive sense defined by the rotation of Y into Z. 

For cylindrical roller bearings, the primary load direction is radial with secondary effects of 
misalignment included. The quasi-static loads applied to the inner ring and shaft and the 
corresponding deflections are represented by four degrees of freedom defined in the inertial 
Cartesian coordinate system. The four displacement degrees of freedom are large Y and Z 
translations and small BY and 0Z rotations. Thus small moments must be generated to enforce 
angular misalignments, but the primary loading is by radial forces. 

The radial displacement of the inner ring and shaft defined in each azimuthal plane is 
Au r = AY cos( 0) + AZsin(0) 

while the corresponding tilt rotation in each azimuthal plane is obtained from a coordinate 
transformation based on the small-angle assumption: 

AB = 0Zcos(0) - BY sin(0) 

This equivalent small tilt angle (A 0) is measured with the positive sense defined by the rotation of 
X into R. Radial translations and tilt rotations of the outer ring and carrier and of each rolling 
element are defined directly in cylindrical coordinates of the inertial reference frame. 

Kinematic constraints leading to rolling element forces applied to inner and outer races are 
developed from large-angle trigonometry in each azimuthal plane of the inertial reference frame. 
Roller/raceway interference is determined in the inertial radial directions leading to radial loads 
only. Roller crown drop is assumed to be in the inertial radial direction although actually defined 
normal to the roller. Maintaining the SHABERTH Hertzian contact equations, while using only 
inertial radial forces and displacements, simplifies the quasi-static analysis. 
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5. BALL BEARING KINEMATIC CONSTRAINT RELATIONS 


Constraint relations for the flexibility enhanced kinematics with angular contact ball bearings will 
be described in this section. The kinematic constraint relations vary only slightly for the case of a 
duplex set. Two similar ball bearings installed back-to-back, as shown in Figure 5-1, are 
kinematically equivalent except for a sign change in the contact angles. For convenience in the 
following discussion, the bearing with the positive contact angles will be called the “primary” 
bearing, while the bearing with negative contact angles will be called the “duplex” bearing. As 
shown in Figure 5-1, the origin of the global inertial coordinate system ( CS t ) is fixed at the initial 
location of the outer race curvature center ( ORCC) of the primary bearing. The ORCC of the 
duplex bearing is displaced along the inertial X axis a distance X cd , which is restricted to positive 
values. Relative elastic deflections of the shaft between the two bearings are represented in the 
structural model of the shaft and inner rings. 



Figure 5-1 Inertial Coordinates for Duplex Bearing Set (Unloaded Condition) 


The basic constraint relations for a ball in contact with both races are first developed without axial 
preload effects. Then modifications of the basic constraint relations to include preload effects for 
a duplex set are discussed. 

5.1 Contact with Both Races 

The following development is applicable to a single ball bearing or to the primary ball bearing of a 
duplex set, only because of an assumed positive initial contact angle. Otherwise the equations are 
equally applicable to either the primary or the duplex ball bearing. 

For each ball location, the kinematic constraint relations are obtained from two different 
expressions for the location of the ORCC. One expression locates the ORCC directly in CS,. The 
other expression locates the ORCC indirectly relative to the location of the inner race curvature 
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other expression locates the ORCC indirectly relative to the location of the inner race curvature 
center ( IRCC) in the inner ring moving coordinate system (CS 2 ) Both expressions represent the 
ORCC radial and axial locations in CSi when the ball is in contact with both races under load. 



Figure 5-2 Mounted Ball-Race Configuration (Unloaded) 


mounting from the ORCC and CS t is seen in Figure 2 to be 


An auxiliary equation may be 
obtained from the two expressions 
of the ORCC radial inertial 
coordinate where the ball is 
unloaded but in contact with both 
races: 

Ro = Rj- BD cos (a 0 ) 

This expression may be rearranged 
as follows for a convenient 
identity: 

BDcos(a°) = Ro-R, (5-1) 

This kinematic relation is shown 
graphically in Figure 5-2. 

The corresponding equation for 
the ORCC axial inertial coordinate 
is null, because both curvature 
centers are initially in the X-R 
plane and the ORCC does not 
move during mounting. The axial 
distance by which the IRCC and 
CS 2 are initially displaced during 


Ae = BD sin(a°) 


(5-2) 


The initial distance between the ORCC and the IRCC in the unloaded mounted configuration is 
seen in Figure 2 to be 

BD = r 0 +r,-D (5-3) 


And the initial (positive) contact angle for the primary bearing is related to diametral clearance by 
a° = cos 1 [(2 BD - P d - AP d )/2BD] (5-4) 

The initial room-temperature diametral clearance is modified by thermal effects to produce a 
revised initial contact angle using an existing analytical method [1], 
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When a ball at the j 1 * 1 azimuthal location is loaded with applied static and centrifugal forces, the 
contact angles at the inner and outer races differ. The j** 1 axial ORCC location under load is the 
axial elastic displacement determined directly in global inertial coordinates as &c OJ . 

The same axial coordinate expressed indirectly in terms of rigid-body and elastic displacements of 
the inner ring and ball results in the following kinematic relation: 

dx oj = Ae+ Ax + 5x,j 

- [r,- (D/2 - 5d tJ )] sin(a, y ) - [r a - (D/2 - &&,)] sin(a (?/ ) 

+ [Oy s\T\(0j) - Oz cos(0j)]{R t + dry 

- [r,- (D/2 - 5dij)\ cos (ay) - [r a - (D/2 - Sd a j)\ cos(a o; )} (5-5) 

The small angle assumption was used with Oy and Oz for transforming quantities defined in CS 2 
into CSi for equation (5-5). Substituting equation (5-2) and rearranging equation (5-5) results in 
the basic axial kinematic relation with the loaded ball in contact with both races: 

0 = BD sin(a°) + 8X } 

- [/*, - (D/2 - sin(a,j) - [r 0 - (D/2 - 5 d 0 j)\ sin(Ooy) (5-6) 

As seen from equations (5-5) and (5-6), some of the rigid-body and elastic displacements are 
collected for convenience into 

SXj = Ax + 8xij - 5x 0 j 

+ [Oy sin(<2^) - Oz cos(<^r)]{/?, + 8r, } 

- [r,- (D/2 - 8d v )] cos (a 0 ) - [r 0 - (D/2 - 8d 0 j)] cos(a q/ )) (5-7) 

The corresponding kinematic relation in the radial direction is determined similarly. The j* radial 
ORCC location under load is the radial elastic displacement (8r OJ ) added to the unloaded location 
( r 0 ). The same inertial radial coordinate expressed indirectly in terms of rigid-body and elastic 
displacements of the inner ring and ball results in the following kinematic relation: 

Ro + Sr OJ = R,+ Ay cos(tf)) + Az sin(<^) + Sr 0 

- — (D/2 - <5t/y)] cos(a !/ ) - [r a - (D/2 - 8d OJ )] cos(a oy ) 

+ [ Oz- cos( 0j) - Oy sin( 0/)] { 8x y 

- [r, - (D/2 - 6 4)] sin(a,y) - [r a - (D/2 - dd OJ )] sm(a OJ ) } (5-8) 

The small angle assumption was again used with Oy and Oz for transforming quantities defined in 
CS 2 into CS, for equation (5-8). Substituting equation (5-1) and rearranging equation (5-8) 
results in the basic radial kinematic relation with the loaded ball in contact with both races: 

0 = BD cos(a°) + SRj 

- [r, - (D/2 - 8d,j)] cos(aj,) - [r 0 - (D/2 - 8d OJ )] cos(a a j) (5-9) 
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As before, some of the rigid-body and elastic displacements are collected from equation (5-8) into 
SRj = Ay cos (0,) + Az sin(<P,) + Sr tj - 5r oj 
+ [Qzcos(<Pj) - 9y sin(<^)] [Sx v 

- [r,- (D/2 - H)] sin(ay) - [r a - (D/2 - 5d 0 J] sin^,)} (5-10) 


The planar kinematic relations given by equations (5-6) and (5-9) are expressed graphically in 
Figure 5-3 as done previously by Davis and Vallance [2], Hams [3] and Jones [4], The latter 
terms in equations (7) and (10), not found in the literature, should be included as they may be 
comparable in magnitude to the ORCC and IRCC elastic displacements which are included. 


IRCC (loaded) 



5.2 Contact with Outer Race Only 

In the unloaded area of its orbit, the f ball can lose contact with the inner race as it is compressed 
against the outer race by centrifugal force. For this condition, the Hertzian contact deflection 
between the ball and inner race (5d t j) is zero. Also for this condition, the outer race contact 
angle (a OJ ) and the outer race Hertzian contact deflection (Sd OJ ) can be determined directly from 
the centrifugal force components in CS 2 . With these restrictions, equations (5-5) and (5-8) are 
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then the limiting relations for lost inner race contact, and the corresponding inequalities for 
defining lost contact in the axial direction is 

5x OJ > Ae+ Ax + Sx,j 

- [r,- ( D/2 - H)1 sin(cXy) - [r 0 - (D/2 - 54,)] sin(a 0 ,) 

+ [Oy s\n(&j) - Oz cos(&j)]{R, + 5r tJ 

- [r,- (D/2 - SdtJ)] cos(a„) - [r a - (D/2 - Sd OJ )] cos(a oy )} (5-1 1) 

and in the radial direction 

R 0 + Sr oj > R,+ Ay cos <<J>,) + Az sin(<Z>/) + Sr v 

- [r,- (D/2 - Sd tJ )] cos (of) - [r 0 - (D/2 - <%,)] cos(a o; ) 

+ [ Oz cos( &f) - 9y sin( &j)] { dx, } 

- [r, - (D/2 - 54)] sin(a y ) - [r a - (D/2 - 54,)] sin(a o; ) } (5-12) 

Substituting equation (5-2) and rearranging inequality (5-11) to isolate terms involving the 
indeterminate inner race contact angle (a,j) results in 

X cj < (r , - D/2) sin(a y ) + [Oy s\n(0j) - Oz cos(<^)] (r,- D/2 ) cos(a y ) (5-13) 

where 

X cj = BD sin(a°) + Ax+ Sx,j - Sx OJ - [r a - (D/2 - 54,)] sin(a o; ) 

+ [Oy sin (®j) - Oz cos (<#)]{/?, + Sr v - [r a - (D/2 - Sd OJ )] cos(a 0 ,)} (5-14) 

Similarly substituting equation (5-1) and rearranging inequality (5-12) to isolate terms involving 
the inner race contact angle results in 

Rcj < (r t - D/2) cos(a y ) + [Oz cos (<ty - Oy sin(<#)] (#v- D/2) sin(a,,) (5-15) 

where 

Rcj = 5Dcos(a°) + Ay cos (4) + Az sin(4) + Sr, } - 5r oj 

- [r 0 - (D/2 - 5d 0 j )] cos(ao,) 

+ [Oz cos (<Pj) - Oy sin(<^')] {&,, - [ r 0 - (D/2 - 54/)] sin(a o; )} (5-16) 

Squaring both sides of inequalities (5-13) and (5-15) and adding the two inequalities eliminates 
the indeterminate inner race contact angle and results in the desired constraint relation: 

X C J + Rcj 2 ^ (r, - D/2) 2 {1 + [Oy s\n(<Pj) - Oz cos(<ty)] 2 } 

When this inequality is satisfied, the j 1 * 1 ball has lost contact with the inner race. 
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The centrifugal force F cj on the f' ball is defined in CS, and directed along the positive r-axis. 
When this force is transformed into CS 2 , using the small angle assumption, the axial and radial 

components are 

F cxj = F C j [Oz cos(<#) - 9y sin(<^)] ( 5 ’ 18 ) 

f.=f (5-19) 

r CTJ r CJ 

For the condition when the ball loses contact with the inner race, the outer race contact angle is 
obtained directly from 

tan(a o; ) = F cx JF cr] = 6z cos (<#) - Oy sin (<#) ( 5 ' 2 °) 

Similarly, when the ball loses contact with the inner race, the Hertzian contact deflection between 
the f ball and the outer race is obtained directly from 

Sd 0 j = (Fcj/Kojf 3 (5 ' 21) 

Because of the small angle assumption used in the coordinate transformations, the magnitude of 
the centrifugal force is slightly larger in CS 2 than it is in CSi. The resulting error is considered 
insignificant. 

5.3 Modifications for Duplex Bearings 

As discussed previously, the primary modification in the constraint relations for the duplex 
bearing is the change in sign of the contact angles. Note that the subscript d is used for variables 
related to the duplex bearing. 

Ordinarily, des-criptive geometrical and mass data for the two bearings in a duplex set are 
identical. Thus the initial (negative) contact angle for the duplex bearing is obviously identical to 
equation (5-4) except for a sign change, and is determined from 

<x° d = -cos'[(2BD d - P dd - AP dd )l2BD d ] ( 5 ‘ 22 ) 

The negative initial contact angle for the duplex bearing in equation (5-22) results in a negative 
value for the mounting axial displacement defined as the axial distance between ORCC and / C 

and given by 

Ae d = BD d sin(a°<) (5 ' 23) 

This equation is, of course, identical to equation (5-2) except that the data relates to the duplex 
bearing. Note that oc^and Ae<j are both negative. 

The only other modification besides contact angle and descriptive data regards the lateral rigid- 
body displacements for the IRCC of the duplex bearing. The five DOFs describing the rigid-body 
translations and rotations of the CS 2 origin are Ax, Ay, Az, Gy and 0z. The axial translation and 
two rotations are identical for the two bearings. However, the IRCC of the duplex bearing is 
displaced in the axial direction relative to the CS 2 origin. 
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(5-24) 

(5-25) 


Therefore, the lateral displacements of the duplex IRCC differ from Ay and Az as follows. 

Ayd - Ay + (X c d - Ae + AeJ) Qz 
Azd = Az- ( X C d - Ae + AeJ) 0y 

With a 0 ,# Ayd, Aid and descriptive data for the duplex bearing, the kinematic constraint relations 
described previously relative to the primary bearing are completely valid for the duplex bearing. 
Based on a°a the loaded contact angles, CLodj and CL,dj, are both negative. 


5.4 Axial Preload Effects 

The extension of the preceding development to account for axial preload effects for a duplex set 
follows the approach of Harris [3] in which all ball forces and contact angles for a bearing are 
assumed equal. Note that the axial components of the preload ball forces in the primary an 
duplex bearings must be equal and opposite, but the radial ball force components and preload 
contact angles could differ. It may also be noted that axial preload effects for single ball bearings 
are accommodated merely by adding the preload force to the applied axial thrust force. For 
duplex sets, an assumption of equal ball forces and contact angles would be valid for the usual 
case where the primary, duplex, and carrier supports are "mirror images" as shown in Figure 5-1 
The two limitations of the present preload analysis are that the number of rolling elements in the 
primary and duplex bearings are equal and that interference fits are not permitted between outer 

ring(s) and carrier. 

The two expressions of ORCC radial location, including elastic deflections due to axial preload, 
produce the following constraint equation for the primary ball bearing: 

R 0 + Sr po = R, + dr pl - (BD + Sn) cos(a^) ( 5 ‘ 26 ) 

Similarly, the two expressions of ORCC axial location after axial preload application produce the 
following constraint equation with elastic structural deflections included: 

Sxpo = BD sin(a°) + Ap + Sx p , - (BD + Sn) sin(a p ) (5-27) 


The ball Hertzian contact deflection (Sn) due to axial preload is related to the specified total 
preload force by the following expression developed by Jones and presented by Harris [3]: 

,o.5 (5-28) 


pre 


n sin(a ) 


= KdD 


BJ 


For convenience, the axial deflection constant (Kd), in units of psi, can be expressed in terms o 
total curvature (5) and elastic constants of the raceway and ball by the following relation where 
Cd is determined from a least-squares fit of the data presented by Harris. 


Kd 


r 

B 
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(5-29) 
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where 


C d = 0.61884472 + 20.2612675 -286.1 5792^ 

+ 2306.1467^-9258.9168#' + 14472.024# 5 


Here B is allowed to vary from 0.0 to 0.2. Note that the elastic constants for both raceways must 
be identical for this representation to be valid. 


Solving equation (5-26) for ball deflection and substituting into equation (5-28) results in the 
following expression involving preload contact angle as the only unknown. 

f - V* 


pre 


»K d B> 


Yl 


= sin(aJ 


5m, 


rp 


cos(a_) 


- BD 


(5-30) 


where, with the identity of equation (5-1), 

5urp = BD cos(a°) +Sr pi - drpo 


(5-31) 


Equation (5-30) may be solved iteratively using the Newton-Raphson procedure for given values 
of elastic structural deflections. In closed form, the equation to be satisfied is 


ctp Op _/( &p)^f ( &p) 


(5-32) 


where 


f(oi p ) 


pre 


6 m , 




nK,B* l D / 


* - sin(a 4^b 


- BD 


/(a ) = -cos(a p i 


5m 




rp 


cos(a ) 


BD 


5m, 




rp 


- ISS^tan^a os(a 


- BD 


J 


A value of a p is determined from equation (32) when a p ' equals a p with sufficient accuracy. For a 
given dp, the corresponding axial and radial ball force components due to axial preload can be 
calculated directly: 


Fxp ~ Fpre/n 
F rp = FVtan(a p ) 


(5-33) 

(5-34) 


The elastic structural deflections of the inner and outer race curvature centers are calculated from 
the appropriate linear elastic flexibility influence coefficient matrices using the ball force 
components given by equations (5-33) and (5-34). Since the fundamental assumption of the axial 
preload calculation is that all balls behave the same under the applied axial preload, the desired 
elastic structural deflections (i.e., 6 V , 5 r ph b Xpo , bx pi ) are defined as the average values of the 
corresponding deflections calculated from the flexibility influence coefficient matrices. T e 
average radial structural deflections are used to update 5 Urp (equation 5-31) for the next iterative 
solution of a p (equation 5-32). When the preload contact angle has been determined with 
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sufficient accuracy, the axial displacement of the inner ring due to preload can be determined from 
equation (5-27) using the expression for ball deflection from equations (5-26) and (5-31). 

Ap = bxpo - 5x pi + 5u rp t?Ln(a p ) - BDsm(a°) (5-35) 

The axial displacement (Ap) of the primary bearing due to axial preload is the essential result of 
the preload calculations described with equations (5-26) through (5-35). The only change in the 
basic kinematic constraint relations is the redefinition of the initial axial distance between CS 2 and 
CSj. Analogous to equation (5-2), this initial axial distance due to mounting and preload is 

Ae p = BD sin(a°) + Ap (5-36) 

The kinematic constraint relations, when the ball is in contact with both races under load, are once 
again obtained from two different expressions for the inertial location of the ORCC. With axial 
preload, the ORCC is, of course, initially displaced axially by the distance bxpo due to preload, 
whereas the IRCC is displaced the additional distance bx pl . However, if the axial elastic 
displacements (bx OJ and bx tJ ) are calculated with ball forces which include the preload force 
components, then the preload elastic structural deflections (6x ,*, and bx pi ) will be included in 5x OJ 
and bxjj. And the gap values defined initially, prior to preload, may be used unchanged with the 
application of axial preload, since changes in the gaps due to preload deflections are determined 
automatically. Similarly, the radial elastic displacements calculated with total ball forces will 
include the preload radial structural deflections (Sr,*, and 5 r p ,) as well. 

Accordingly, the only change in the previously developed kinematic constraint relations, for the 
primary bearing of a duplex set, is to replace Ae (equation 2) with Ae p (equation 5-36). Equation 
(5-6) remains unchanged while equation (5-7) has only one added term (Ap) to account for the 
preload axial displacement: 

SXj = Ax + Ap + dxjj- Sx 0 j 

+ [Gy sin (&j) - Gz cos( <£?)]{./?, + Sr v 

- [r,- (D/2 - tidy)] cos(ajj) - [r a - (D/2 - Sd Q j)] cos(a 0 ,)} (5-36) 


The radial kinematic relations given by equations (5-9) and (5-10) are unchanged. 

For a primary bearing, the analogous discussion applies to the constraint relation for a ball 
contacting the outer race only (inequality 5-17). For the axial component, Ae is replaced by Ae p 
so that equation (5-14) becomes 

X C j = BD sin(a°) + Ax + Ap + Sxy - Sx OJ - [r a - (D/2 - 5 d OJ )] sin(a o; ) 

+ [Gy sm(<Pj) - Gz cos (<#)] {R, + Sr v - [r a - (D/2 - Sd OJ )} cos(a 0 ,)} (5-37) 

Equation (5-16) for the radial component and inequality (5-17) are unchanged. When inequality 
(5-17) is satisfied, the j* ball has lost contact with the inner race. 

Similar modifications are required to incorporate axial preload effects for the duplex bearing. 
With negative contact angles and with reversed axial ball force components for the duplex 
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bearing, the magnitude of the axial displacement of the inner ring due to axial preload becomes 
analogous to equation (35): 

Ap d = 5x pdo - bxpd, + &u rpd tan(a pd ) - BDs\n(a% ( 5 * 38 ) 

Here 6 *,„ and a° d are obtained from equations (5-28) through (5-32) using descriptive 
geometrical and mass data for the duplex ball bearing set. Ordinarily, both bearings of a duplex 
set will be identical, but the only limitations in the preload calculation, as noted, are that both 
bearings have the same number of rolling elements and that outer nng(s) and earner have no 

interference fit. 

Analogous to equation (5-23), the initial axial distance between CS 2d and CS, increases by the 
quantity 

Ae dp = BD d sin(a° £ ) + Ap d (5 ' 39) 

Of course, for the duplex ball bearing, Ae dp (equation 5-39) must replace Ae in equations (5-2) 
and (5-14) to give kinematic relations for the duplex bearing analogous to equations (5-36) and 
(5-37) for the primary bearing. 

Also, for the duplex bearing, equations (5-24) and (5-25) must be changed to redefine the lateral 
displacements of the duplex IRCC ( Ay and Az) as follows: 

Ay dp = Ay + (X ed - Ae p + Ae dp ) 9z (5 ' 40) 

Az dp = Az - {X cd -Ae p + Ae dp ) 6y (5 ' 41) 

These kinematic constraint relations are consistent with the prescribed axial preload force 
Therefore, ball forces resulting from these constraint relations will properly include the 
corresponding preload force components. The radial preload force components are se - 
equilibrating as are axial preload force components for duplex bearings. For single ball bearings, 
as noted previously, the axial preload force is merely added to the externally applie axia orce 

component. 

5.4.1 Requirements 

Use of the preload spring option in FEREBA must conform to the following requirements: 

1 . Axial preload springs may be selected for most ball bearing configurations (BSCOR = 
2,3,. ..,6) to represent distributed macroscopic preload effects. For the single ball 
bearing, an axial preload spring is installed between the outer ring and carrier on one 
side of the outer ring. For the duplex ball bearing set, an axial preload spring is 
installed between the two ball bearings. For both configurations, the axial preload 
force on the outer ring is transmitted through the rolling elements to the inner ring and 
shaft as shown in the attached sketch. Note that there is positive clearance between 
the outer ring and carrier for the bearings being considered. An interference fit is 
permitted only between inner ring and shaft, not between outer ring and carrier when 
axial 'preload is involved. 
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2. The following local loading capabilities are required for the idealized preload spring: 

a. completely unloaded locally such that no force is developed; 

b. compressed locally within a linear range such that a compressive force is 
developed; 

c. completely bottomed locally such that the full compressive constraint force is 
developed. 

3. The main function of the preload spring is to impose an axial force and two moments 
on the outer ring. These applied macroscopic preloads remove all bearing free play 
and increase the bearing stiffness. In addition to these effects of discrete axial forces 
applied to the outer ring, corresponding Coulomb friction forces are developed which 
minimize outer ring tangential rotation. 

5.4.2 Modeling Approach 

The modeling approach is the simplest preload spring model which meets the above requirements, 
uncoupled axial linear springs at each ball location. This model provides the three required 
loading capabilities (unloaded, loaded and bottomed) as well as the macroscopic preload total 
force and moments. The total preload spring force is simply the sum of the individual forces, and 
the total moments are simply the products of the individual forces times their respective moment 
arms. An individual force could be zero, the linear spring force, or the appropriate constraint 
force, depending on local (and global) kinematic conditions. 

This model does not represent theoretical stiffness coupling occurring around the circumference in 
some types of actual preload springs, nor does it represent nonlinear behavior resulting from 
Belleville springs. Since actual preload springs under consideration have stiffness terms on the 
order of one percent of the outer ring stiffnesses, any error in the outer ring deflection due to 
stiffness coupling or nonlinear behavior would be negligible. Local deflections of the preload 
spring will be defined essentially by the axial displacement and tilt of the outer ring relative to the 
carrier or to the other outer ring. The relatively soft preload spring linear stiffnesses will have a 
small effect on outer ring deflection so long as bottoming of the preload spring is properly 
represented. 

The basic assumptions used in representing the preload spring forces are as follows: 

a. The continuous preload spring is divided into n independent axial springs contacting 
the outer ring at the n ball azimuthal locations. These separate springs are combined 
with the finite-element model of the outer ring and carrier generated by ANSYS. This 
combination of preload spring stiffnesses with the ANSYS stiffness matrix is 
accomplished in FEREBA prior to the matrix condensation and inversion operations. 

b. The prescribed total preload force, which is user specified, is developed by the 
individual axial preload springs being compressed equally in their linear range. 

c. The configuration corresponding to the mounted, preloaded ball bearing system has 
inner ring, shaft and spacer all rigidly connected with no gaps separating these 
components. This configuration results from an interference fit of the inner ring on the 
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shaft and from the axial locking of the inner ring and spacer against the shaft shoulder 
during preloading. 

d. All initial clearances between the outer ring and carrier, including the gaps at the 
preload spring locations, are defined before the preload is applied. By definition, the 
initial gap separation distance at a preload spring is the minimum distance consistent 
with zero initial force. Axial displacements required to generate the specified preload 
are used to define the net axial clearances after mounting. Ball forces and bearing 
outer ring radial structural flexibility are used to define the net radial clearances after 
mounting. For single ball bearings, the specified initial preload force is simply added 
to the applied axial thrust force. For duplex ball bearing sets, a preliminary preload 
analysis is performed. This preliminary preload analysis determines ball forces, elastic 
deflections of the inner and outer rings, and initial mounted contact angle; but the basic 
result is simply the axial displacement of the inner ring relative to the outer ring due to 
initial axial preload. 

After mounting, the force in the preload spring is allowed to vary in response to the externally 
applied forces and moments. The specified constant preload force is distributed to the n ball 
azimuthal locations and applied as equal-and-opposite external forces. Preload force changes in 
the linear range due to structural displacements are obtained with the basic solution procedure, 
since the individual axial preload springs are included in the AN5YS model. The basic solution 
procedure also determines the constraint forces when the preload spring bottoms locally. 
Additional logic is required to monitor local lift-off of the preload spring and to cancel the small 
incremental tensile force incorrectly generated by the ANSYS linear model. 



6. ROLLER BEARING KINEMATIC CONSTRAINT RELATIONS 


The basic assumption of the quasi-static roller bearing analysis is that only forces in the inertial 
radial direction are applied by each roller to the inner and outer raceways. Thus, although small 
moments may be applied to the inner ring and shaft corresponding to angular (tilt) misalignment, 
the primary loads are radial, and axial thrust loads are not permitted. 

The classical roller bearing analysis by Jones [2] permits only radial loads and displacements with 
no angular misalignment. Subsequent developments by SKF, as published by Harris [3] and Liu 
[6], include the effects of angular misalignment (tilt) in approaches similar to that employed by the 
SHABERTH computer code [7], More comprehensive analyses, such as that recently developed 
by SKF [8], include the axial degree of freedom for all types of roller bearings in a formulation 
similar to that previously used for ball bearings. However, the scope of the present frictionless 
quasi-static analysis for cylindrical roller bearings is limited to the radial and tilt freedoms 
consistent with SHABERTH, which serves as the development platform for this effort. 




A definition of basic geometrical relations is provided in Figure 6-1 for the unmounted roller 
bearing configuration. As shown, the radial distance from the unmounted bearing center to the 
outer raceway is 

R < , = DJ2 + D/2 + Pj/2 ( 6_1 ) 

The radial distance from the bearing center to the inner raceway is 

R, = Rt-D-PJl < 6 * 2 ) 

As used in SHABERTH, the analytical expressions for roller/raceway contact forces are defined 
in terms of kinematic interference between roller and raceway. These kinematic constraint 
relations are defined in terms of inertial radial displacements and inertial tilt rotations of the roller 
centroid at each azimuthal roller position. Due to the shaft/bearing system modeling approach 
used by SHABERTH, axial displacements of the roller centroid in a rigid body sense were 
permitted. In the present development, the analysis of a single roller bearing eliminates the need 
to consider rigid-body axial motion. Consistent with precluding axial, roller displacements, the 
raceway widths are assumed large compared with the roller lengths. 

6.1 Interference Between Roller and Outer Raceway 

The initial mounted (unloaded) position of roller and outer raceway at the j* 11 azimuthal plane has 
contact between roller and raceway such that points A 0 and B 0 , shown in Figure 6-2 for the 
loaded condition, are initially coincident. Under applied load, the outer raceway in the j azi- 
muthal plane translates in the radial direction a distance bu OJ and rotates in the tilt direction 
through a small angle Go,. 

Simultaneously, the roller 
centroid translates radially 
a distance Ary, and rotates 
through a small angle 
For each roller. Ary = X(2) 
and $ = X(3) are the 
independent variables in 
the FEREBA equilibrium 
solution. The radial inter- 
ference between point A 0 
on the roller and the 
radially adjacent point C 0 
on the outer raceway is 
Sr,,, shown negative for 
clarity in Figure 6-2 as a 
separation rather than an 
interference. 
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With radial translation but no rotation, the radial separation at point A 0 is the distance dj in 
Figure 6-2 given by 

di = Su 0 j-Arj 

The radial increments associated with the two inertial rotations are shown in Figure 6-2 as d 2 and 
ds and are given analytically by 

d 2 = DI 2 sin(^) tan(&/) 
di = D/2[\- cos(^)] 

The total radial interference between points A 0 and C 0 is then 

Sr OJ = d 2 - d, - d 3 = Arj - Su OJ - DI 2[1 - cos($ - e oj )/cos( 0 o j)] (6-3) 

The radial interference at the center of the k* slice of the j* roller (Sqokj) is shown negative in 
Figure 6-2 as a separation rather than an interference for clarity. An analytical expression for t s 
interference, which includes a slight approximation for crown drop, is 

Sq^ = Sr OJ + Xok sin($ - 0 o ^lcos(0 OJ ) - CD k ( 6_4 ) 

The crown drop approximation for the incremental decrease in radial interference is 

CD k * CD k cos($- 0 O j)lcos(0 O j) (6 ' 5) 

This simplifying approximation is, for the small angles expected, accurate to several significant 
digits. Probably the crown drop data, which can include raceway crown in addition to roller 
crown, has significantly greater variability than the error due to this approximation 

The basic kinematic constraint relations for radial interference, equations (6-3) and (6-4), are 
subsequently used to determine the Hertzian contact forces between the roller and outer nng. 


6.2 Interference Between Roller and Inner Raceway 

The initial mounted (unloaded) position of roller and inner raceway at the f azimuthal plane has 
maximum radial clearance between roller and raceway at points^ and B it shown in Figure 6-3 tor 
the loaded condition. Under applied load, the roller centroid in the j azimuthal plane translates 
radially a distance Ar h and rotates through a small angle £. As with the outer nng, Ar p - X(2) 
and <| >, = X(3) are the independent variables in the FEREBA equilibrium solution. 

Simultaneously, the inner raceway translates in the radial direction a distance Aw* and rotates in 
the tilt direction through a small angle AG,. These two displacements of the inner raceway 
combine the imposed inertial translations and rotations of the inner nng and shaft with elastic 
distortions of the inner raceway in the j* azimuthal plane 

Au n = AY cos (0j) + AZ sin(d>y) + Su u ( 6 ’ 6 ) 

A0j = 0Z cos(d> y ) - 0Y sin {(Dj) + 0,j (6 ' 7) 
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The radial interference 
between point Ai on the 
roller and the radially 
adjacent point C, on the 
inner raceway is Sr, h 
shown negative in Figure 
6-3 as a separation rather 
than an interference for 
clarity. With radial trans- 
lation but no rotation, the 
radial separation at point 
Ai is the distance <4 in 
Figure 6-3 given by 

d4 = ATj - AUrj 

+ Pd/2+AP d 

The radial increments 
associated with the two 
inertial rotations are 
shown in Figure 6-3 as d$ 
and da and are given 
analytically by 

di = D/2 (1 - cos($)) 

d 6 = Rr[\ - cos( Qj)\l cos( &j) + DU sin($) tan(0,) 

The total radial interference between points Aj and C, is then 
Sr t = d 6 -d 4 -d } = Aurj-Arj-Pd/2-APj 

+ R,[ 1 - cos( 6j)\l cos( Oj) - D/2 [l - cos(4 - ^)/cos(^)] (6-8) 

The radial interference at the center of the k* slice of the j A roller (£ 7 *) is shown negative in 
Figure 6-3 as a separation rather than an interference. An analytical expression for this 
interference, which includes basically the same approximation for crown drop as used for the 
outer raceway, is 

5qog = Sr v + x* sin(^ - $)/cos (0j) - CD k (6* 9 ) 

The basic kinematic constraint relations for radial interference, equations (6-3) and (6-4) for the 
outer raceway and equations (6-8) and (6-9) for the inner raceway, are subsequently used to 
determine the roller Hertzian contact forces. The equations for roller forces and moments, 
required to calculate elastic deflections of the outer and inner raceways, are developed in the 
following section. 



Figure 6-3 Interference Between Roller and Inner Raceway (Loaded) 
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6.3 Roller Forces and Elastic Displacements 

The roller/outer race Hertzian line contact force on the k* slice of the j* roller is 

Q* K 0 j (Sqokj ) 10/9 


and the corresponding force at roller/inner race contact is 


Qng = K,j (Sq,ig ) 10/9 


(6-11) 


Each slice along each roller is checked independently for interference with both raceways. With 
positive interference, the Hertzian contact force is positive in compression. With negative 
interference (separation), the Hertzian contact force is zero. Qok? rnay be defined as t c 1 
force applied to the outer raceway by the roller. Q ikj may be defined as the radial force applied to 

the roller by the inner raceway. 


' The total radial forces in the f azimuthal plane are obtained by summing the contact forces at 


each slice: 


Poj ~ E Qokj 
P ij ~ E Qikj 


(6-12) 

(6-13) 


Here, again, P OJ is positive for a compressive force applied radially outward by the j roller to the 
outer raceway, while P i} is positive for a compressive force applied radially outward by the inner 
. raceway to the j lh roller. 


The total tilt moments applied to the f roller by the two raceways are similarly obtained by 
summing the contributions of the contact forces at each slice. 

Moj^-IQokjXok (6_14) 

My = ZQikjXtk (6 ' 15) 


These tilt moments are defined positively, by the rotation of X into R, about the roller centroid 


Radial forces and tilt moments applied to the outer and inner raceways are converted to two 
statically equivalent radial forces on each raceway. These two forces are separated by a constant 
distance defined as a percentage of the roller effective length. The two forces applied to the outer 
raceway in the j* azimuthal plane for the right and left nodes of the outer ring finite-element 

model are given as follows: 

F orj = P 0 J2-M 0 J{^L 0 ) (6 ' 16) 


Similarly, the two forces applied to the inner raceway are given by 

Firj — — P ij 2 — A^/(A, Li) 

Ftij = -Pi,- Fa 


(6-18) 

(6-19) 
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The two dimensionless constants Xo and X, in equations (6-16) and (6-18) are set to values of 
0.60, making the length of the moment arm (A, Z,)in equation (6-18) equal to 60 percent of the 
effective roller length. This value was chosen such that the maximum and minimum bending 
moments are approximately equal for a uniformly loaded beam on two supports. The finite 
element model master raceway nodes must be defined at locations established by this distance. 

The forces applied to the inner and outer raceways (equations 6-16 through 6-19) produce 
corresponding elastic deflections at the same locations when applied to the appropriate finite- 
element model. These elastic deflections are then converted to the quantities required for the 
kinematic constraint relations (equations 6-3 through 6-9) as follows: 


du oj = (Sr on + Sr 0 , } )/2 C 6 ' 20 ) 

e oJ = - Sroij)/(Xo L 0 ) ( 6 - 21 ) 

du tJ = {Sr ir] + Sr ,ij)/2 ( 6 ’ 22 ) 

6 t j = (8r ir j - 8ruj)/(Xi L,) ( 6 - 23 ) 


This complete set of equations results in a consistent system of forces and compatible 
displacements. The following sections outline the implementation of these equations into the 
FEREBA code. 
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7. NONLINEAR ELASTIC DEFLECTIONS OF OUTER RING AND CARRIER 


The general finite-element model of the outer ring and carrier for single or duplex ball bearings 
contains soft springs to accommodate deadband and preload springs to accommodate axial 
preload effects. The corresponding flexibility matrix is partitioned as follows after the ANSYS 
stiffness matrix is transformed from ANSYS cylindrical coordinates to FEREBA cylindrical 
coordinates and reduced to eliminate the unnecessary DOFs: 


r ^ 

Uj 


~CO,j 

co 12 

8 

u. 


0 

'U 2 

> — 

C0 21 

co 22 

CO„ 

< 


Vs\ 


C0 31 

CO S2 



.Fc. 


where {Ui} = elastic axial (X) and radial (. R ) deflections of ORCCs, or for a roller bearing, the 
raceway contact points, . 

{t/?} = elastic X and R deflections of rolling element/outer race contact points where 
rolling element forces are applied, 

{ U 3 } = elastic X or R deflections of candidate contact points or gaps between the outer 
ring(s) and carrier, 

{Fbo} - X and R rolling element contact forces applied to the outer race, and 

{F c } = X or R compressive constraint forces applied to contact points between the outer 
ring(s) and carrier when the gaps are closed. 


Note that vectors {£/;}, {U 2 } and {F* 0 } have nodes in specific order for each of the outer 
ring/carrier support configurations defined by the index IBSCOR. Each of the nodes specified at 
the ORCC and ball/race contact points has an X and R DOF. The vectors { U3} and {F c } also 
have specific nodes for a given configuration, but each node has either an X DOF or an R DOF 
but not both. Note that the ball forces {Fjo} and the flexibility coefficients (CO,y) are known, but 
all the elastic deflections and the me constraint forces {F e } are unknown. 


Equation (7-1) may be written as follows with the unused partitions omitted: 


Uj 

U 3 


CO I2 

C0 32 



(7-2) 


In this form, the flexibility matrix is directly applicable to the three roller bearing support 
configurations as well as to the six ball bearing configurations. The variables in equation (7-2) are 
now defined more generally as follows: 

{ Ui } = elastic bearing deflections required for subsequent kinematic relations 

{ U 3 } = elastic deflections of candidate contact points or gaps between the outer ring and 
carrier 

(F*o} = rolling element applied force vector 

{F c } = compressive constraint forces applied to contact points between the outer ring and 
carrier when the gaps are closed 
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[< CO 12 ] = flexibility coefficients relating desired bearing deflections { U/} to rolling element 
forces {Fio} 

\COii\ — flexibility coefficients relating desired bearing deflections {£//} to constraint forces 
{Fc} at closed gaps 

[CO 32 ] = flexibility coefficients relating gap deflections { U 3 } to rolling element forces {F^} 
[COjj] = flexibility coefficients relating gap deflections {U 3 } to constraint forces {F -} 

It may be noted that the partitions of equation (7-2) comprise the entire flexibility matrix for the 
roller bearing configurations. This is because the desired bearing deflections {Uj} are the same as 
the deflections {U 2 } of the points at which the rolling element forces are applied. 

Each gap closure involves only a single axial or radial force magnitude which is normal to the 
contact surface. Friction forces in the plane of the contact surface are neglected, and the DOFs 
• corresponding to friction forces have previously been reduced out of the flexibility matrix. 

The constant contribution to {C/3} is the displacement vector { U 3 b } corresponding to the known 
rolling element forces 

{U 3b } = [C0 3 2]{Fbo} (7-3) 

where matrix [CO32] is of order me by mbo. This constant deflection vector may be calculated 
immediately. The basic linear elastic equation relating constraint forces to the corresponding 
. elastic deflections is 

{U3}-{U 3b } = [COssUFc} (7-4) 

where matrix [C0 33 ] is of order me by me. From the basic matrix relation of equation (7-4), the 
row corresponding to the point on the outer ring at the k* gap is 

U 3 (O k ) - U 3b (O k ) = lc0 33 (0 kj ) * F C (J) (7-5) 

where O k denotes the X or R elastic deflection of the point on the outer ring at the k* gap, andy is 
the index of all me candidate constraint forces. Similarly, the row corresponding to the adjacent 
point on the carrier at the k* gap is 

U 3 (C„ ) - U u (C t ) - T.CO sl (C kJ ) * F C (J ) (7-6) 

where C k denotes the X ox R elastic deflection of the point on the carrier at the k* gap. 

Positive faces are defined as those faces on the outer ring for which positive deflections of the 
outer ring initiate contact with a stationary carrier. By definition of the cylindrical coordinates, all 
radial faces of the outer ring are positive faces. Roller bearings are defined to have only radial 
DOFs. Ball bearings have both radial and axial DOFs. Axial faces on the positive X, side of the 
ball bearing ORCC are also positive faces. Thus, for IBSCOR = 2 in section 3.1, the gaps 
numbered 2, 3 and 4 are on positive faces. Axial faces on the negative Xj side of the ball bearing 
ORCC are negative faces. In the section 3.1 figure, gap number 1 is on the negative face. For 
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duplex ball bearing sets, the face between the primary and duplex bearing is defined as positive or 
negative with the same criterion as for single bearings. In all cases, positive contact forces are 
compressive. 

7.1 Kinematic Constraints for Positive Faces 

When the k* gap of a positive face is closed, the constraint forces on both sides of the gap are 
equal and opposite: 

FAO„)--F,(C l ) ( 7 ‘ 7 ) 

where F c (C k ) is the positive (independent) constraint force applied by the outer ring to the earner 
at the k* gap, and F c (O k ) is the negative constraint force applied by the carrier to the outer ring at 
the k** 1 gap. 

When the k* gap is closed, the points on both sides of the contact surface displace equally. This 
results in the following kinematic constraint relation based on the small deflection approximation: 

U 3 (O k ) = U 3 (C k ) + G k ( 7 - 8 ) 

where G k is the input initial X or R clearance at the k 01 candidate contact point between outer ring 
and carrier. By definition, a positive value of G k indicates clearance, whereas a negative value of 
G k indicates interference. 

For use in testing whether the k* gap is closed, equation (7-8) may be written as 

Ag(k) = Ui(C k ) + G k - U 3 (O k ) ( ? -9) 

Then if Agj k) <0, gap closure may be assumed and the corresponding constraint forces must be 
calculated to enforce Ag(k) = 0. 

A relation corresponding to equation (7-9) in terms of { f/?*} is, for gaps at positive faces, 

Agb(k) = U 3 b(C k ) + G k - U3b(O k ) (7-10) 

The variables Ag b (k) may be calculated immediately since the terms are dependent only on known 
quantities. Subtracting equation (7-6) from equation (7-5) and substituting the relations in 
equations (7-8) and (7-10) gives the constraint equations in efficient computational form for 
positive faces: 

me 

Ag b (k) = £ [C0 33 (0 kj ) - C0 33 (C kj )] * F C (J) (7-11) 

J-J 
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7.2 Kinematic Constraints for Negative Faces 

When the k* gap of a negative face is closed, the constraint forces on both sides of the gap are 
equal and opposite but with signs reversed from equation (7-7): 

F c (C k ) = -F c (O k ) ( 7 - 12 > 

where F c (O k ) is now the positive (independent) constraint force applied by the carrier to the outer 
ring at the k* gap, and F c (C k ) is now the negative constraint force applied by the outer ring to the 
carrier at the k** 1 gap. 

When the k* gap is closed, the points on both sides of the contact surface displace equally. This 
results in the following kinematic constraint relation based on the small deflection approximation: 

U 3 (C k ) = U 3 {O k ) + G k ( 7 - 13 ) 

For use in testing whether the k* gap is closed, equation (7-13) may be written as 

Ag(k) = U 3 (O k ) + G k - U 3 (C k ) ( 7_l4 ) 

Then if Ag(K) < 0, gap closure may be assumed and the corresponding constraint forces must be 
calculated to enforce Ag(k) = 0. 

The relation corresponding to equation (7-14) in terms of {U 3 b} is, for gaps at negative faces, 

Agb(k) - U 3 b(O k ) + G k - U 3 b(C k ) ( 7_ I 5 ) 

Subtracting equation (7-5) from equation (7-6) and substituting the relations in equations (7-13) 
and (7-15) gives the constraint equations in efficient computational form for negative faces: 

4g,(*) - 2 [CO„(C t ,) - CO n (O k j)] * F'U) ( 7 - 16 ) 

j‘I 

7.3 Solution Procedure for Assumed Gap Closures 

Equations (7-1 1) and (7-16) may be written in matrix form as follows for only those mcgc gaps 
assumed to be closed 

{Agb} = 1-T,]*[C0 33 ]*[T,] t *{F cc } ( 7 -17) 

This format relies on the fact that matrix [C0 33 ] is symmetric. Matrix [T s ] is of order mcgc by me 
and is defined to extract vector {Agb } from vector {U 3 b} using equations (7-10) and (7-15). For 
positive faces, matrix [ T s ] is a sparse matrix wherein the k** 1 row consists of all zeros except for a — 
1.0 in the column corresponding to point O k and a +1.0 in the column corresponding to point C k . 
For negative faces, the signs of the coefficients corresponding to points O k and C k are reversed. 
Note that the sign of matrix [Y,] must be reversed, as shown, for extracting the proper influence 
coefficients using equations (7-11) and (7-16). 

The modified force vector {Fee} contains only mcgc non-zero independent constraint forces 
defined in equations (7-7) and (7-12). . Note that there is only one independent constraint force 
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per gap and that this constraint force is positive in compression. The relationship between {F C c} 
and {Fc} is expressed in matrix form as 

{F c } = |X] T *{F«} (7-18) 

The independent constraint forces required to satisfy the kinematic constraints of equations (7-8) 
and (7-13) for the mcgc gaps assumed to be closed are then calculated simply by inverting 
equation (7-17) 

{Fee} = {[-Ts]*[COssnTsfy l *{Ag b } (7-19) 

The inverted matrix is generally not positive definite, but it is non-singular. The validity of the set 
of assumed gap closures could be assessed by checking for impossible negative (tensile) in- 
dependent constraint forces {F cc } and impossible negative (overlapping) gap separation distances 
{4g}. Equation (7-20) for {Ag} is the matrix formulation of equations (7-9) and (7-14): 

{Ag} = [Tf]aU 3b } + [COv]{F c }) + {G} (7-20) 

where [7/] is the full matrix [7’,]. The check relation actually used is slightly more complex than 
equation (7-20) because of the nonlinear preload spring which has zero force when unloaded. 


7.4 Preload Spring Unloading 

For single ball bearings or duplex ball bearing sets which admit axial preload, the preload spring is 
represented by appropriate axial stiffness terms at each azimuthal plane. These local preload 
springs are treated like the fictious soft springs required to convert the stiffness matrix of the 
outer ring(s) and carrier into a flexibility influence coefficient matrix. Therefore, the basic analysis 
procedure will properly represent the preload force in each local spring both in the linear spring 
range and after the preload spring has bottomed. However, when the separation distance at one 
of the preload spring gaps exceeds the initial G* separation, the basic analysis procedure will 
calculate an impossible tensile force in this local preload spring. Some modifications to the basic 
procedure are then required to null the preload tensile force. This modification may be 
summarized as follows for the conditions relating to a local axial preload spring. 

At each azimuthal plane, the local preload spring is defined to have zero force when the initial 
spring length equals G k . When the separation distance at a preload gap decreases from its initial 
value, a compressive force is properly developed in the preload spring. When the separation 
distance at a preload gap increases from its initial value, an impossible tensile force is developed 
which must be nulled by a self-equilibrating pair of externally applied constraint forces with 
magnitude F pre ,. When the relative deflection between the two surfaces of a preload gap exceeds 
zero Ag p > 0, the required compressive constraint force required to null the impossible tensile 
force is 

F pre , = Kpre, * Agp ( 7 ‘ 21 ) 

where K pre , is the linear stiffness coefficient of each local preload spring. For each ball bearing 
configuration which accepts axial preload, the appropriate signs are given to F pre , to form the pair 
of external constraint forces contained in vector {F cp } . These additional constraint forces are then 
included in the solution procedure along with those constraint forces due to closed gaps obtained 
from equation (7-18). 
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8. PROGRAMMING IMPLEMENTATION 

This section is provided as a guide to installing and maintaining the FEREBA software. Sections 
are also included that describes the nonlinear iteration methodology incorporated in the program. 
For reference, a flowchart is provided which lists all the subroutines in the code. 


8.1 Flowchart 

A flowchart displaying the subroutine names and execution hierarchy is contained in Figures 8-1 
and 8-2. In the first figure, routines that read the analysis file (READAF module), perform 
problem initialization (CALCON module), and create the compliance matrices (FEMFLX module) 
are listed. 


FEREBA 

|— DBLSE T sets double precision arrays to a constant (used in many routines) 

| — READAF reads analysis variables and bearing data from analysis file written by PREBAN 

T— GETARG returns command line arguments (note - system specific routine) 


I— DATOUT 

t FDATE 
ROLDAT 


prints analysis summary - geometry, loads, fits, materials, temperatures 
, writes analysis time and date to output file (note - system specific routine) 
prints data on rolling elements specific to ball or roller bearing 


| — CALCON calculates general constants relevant to all bearing types 

|— BCON defines parameters specific to ball bearings 

f— ABDE L calculates dimensions and elastic constants for point contact 

I— CRCON defines parameters specific to roller bearings 

I— SLI CE S calculates slice constants and crown drops for cylindrical roller bearing 

— SE TGAP sets gap between roller flange and roller end depending on flange type 

I — ABDEL calculates dimensions and elastic constants for line contact 

I — SPRING determines approximate spring compliances for axial, radial, and angular loads 
— INDEL makes initial guess on inner ring global deflections 
I — TMPFI T calculates change in clearances due to differential temperature effects 


— RD8ASC 
— MATPRT 
— TRANAF 
— ADDSTF 
—HATRED 


| — FEMFLX reads finite element substructures and converts to reduced flexibility matrices 

-ETIME elapsed time counter (note - system specific routine) 

reads ASCII format ANSYS substructure matrices, stonng stiffness only 
prints out matrices 12 columns per page (for diagnostics) 
reorders DOFs in stiffness matrices to FEREBA system 
adds 'dummy' springs for preload & deadband 

condenses DOF out of cylindrical stiffness matrices 

lllsH I FTS shifts rows and columns of stiffness matrices over removed DOF 
— DLFI RG IMSL routine to obtain solution to linear equation set 
— DLFTRG IMSL routine that computes LU factor of a matrix 
—COD IM sets dimensions of outer ring/carrier flexibility matrix 
— COART stores required partitions of outer ring/carrier flexibility matrices 
— CIDIM sets dimensions of inner ring/shaft flexibility matrix 
— C I PART stores required partitions of inner ring/shaft flexibility matrices 


continues 


Figure 8-1 FEREBA Flowchart (Initialization & Matrix Manipulation Routines) 


[FEREBA] 

|— QS TATE drives solution of quasi-static equilibrium problem 
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-ETIME elapsed time counter (note - system specific routine) 

-GAPDE F sets values of gaps depending on ring configurations 

-PRELOD determines axial displacement of inner ring due to axial preload 

l— NRSOLV determines axial preload contact angle using a Newton-Raphson iteration 

l— DE FLOP calculates average linear elastic ORCC deflections due to axial preload 

I — DE FLIP calculates average linear elastic IRCC deflections due to axial preload 

-DO 1 , NITER iteration loop on elastic displacements (NITER set by user) 

SLVBEQ directs solution of bearing inner ring force equilibrium equations 

—DO NITRF iteration loop on inner ring force equilibrium (NITRF set by user) 

If— DO 2, NBRG loop on each bearing (maximum of 2) to get loads or deflections 
L- GUE SS get initial guesses of cage and ball/roller unknowns 

GROLLB sets roller positions, rotational and orbital speeds 

GBALLB sets ball positions, rotational and orbital speeds 

I GUE SCG estimates cage-ring relative motion & assigns cage DOFs* 

— BRGFOR determines bearing loads to match assumed or given deflections 
[— PRE PAR initialize rolling element, cage variables, and arrays 
1— BALLIN calculates curvature center distances (ball bearing) 

— UNLODB determines if ball is out of contact with inner race 

— UNLODR determines if roller is out of contact with inner race 
I — MAXMIN sets max/min values for cage and rolling element variables 

— DNEQNF IMSL solution for nonlinear equations with numerical partials 

>— BE ARE Q sets up equilibrium equations depending on bearing type 

1— BALLEQ defines ball bearing equilibrium equations 

t FMIX calculates frictional load components* 

BDRAG calculates drag coefficient* 

— RQLLIN transforms coordinates for interference calculations 

T— FLNDE F calculates roller/flange interference 
I — ROLLEQ defines roller bearing equilibrium eqns 
— TPNORM calculates loads about roller CG 

— FLNORM determines flange loads 

— FMIXR calculates frictional load components* 

I — FLMIX calculates flange frictional loads* 

— BRGAX calculates rolling element acceleration terms* 

I — CAGEEQ defines cage-element interaction* 

— CALFOR calculates forces acting on outer ring 

' — SUMFI R sums forces and moments acting on inner ring 

— DLST.RG IMSL routine to solve set of linear force equilibrium equations 

I — GUE SS get final values of cage and ball/roller deflections and speeds 

— DEFLOR calculates elastic deflections of outer ring/carrier 

I— CNCLSD produces a set of assumed gap closures (based on current deflections) 

— CNSTI F computes constraint transformation matrix (for gap closures) 

— LFTOFF modifies constraint forces to null tension 

— DLINRG IMSL routine to invert constrained flexibility matrix 

— DEFLIR calculates elastic deflections of inner ring/shaft 
— e RRCAL checks values of elastic deflections, revises if necessary 
—LI FE calculates bearing fatigue life 

1— FLMFAC lubrication film factor for life adjustment* 


I — OUTRSP writes solution data to print file 
f— RE LOUT writes rolling element output 

I — RE FLOT writes output for roller bearing flanges — 

Figure 8-2 FEREB A Flowchart (Main Calculation Section Routines) 

In Figure 8-2, which illustrates the primary quasi-static equilibrium modules, several routines are 
marked with an asterisk (*). These routines are not required for the current version of FEREB A, 
and were included for future editions which intended to add friction to the calculations. In the 
present code, these routines are simple shells, and perform no functions. 
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8.2 Nonlinear Elastic Deflection Implementation 


The nonlinear elastic deflection solution is one of the more significant features of FEREBA. e 
solution procedure described in section 7 is implemented in subroutine DEFLOR usmgro mg 
element forces {F^} and initial gap separation distances {G} input from subroutine QSTATE and 
using flexibility influence coefficient partitions. The output of DEFLOR is the elastic deflections 
{U lb }, designated {Gy} in equations (7-1) and (7-2), of the ORCC or roller bearing desired or 
subsequent use in the kinematic constraint relations. Subroutine DEFLOR is called iteratively 
until the rolling element forces {F^} and the corresponding elastic deflections {U lb } are 
consistent with the kinematic constraint relations. 


The solution procedure leading to outer ring elastic deflections with deadband is accomplished in 
subroutine DEFLOR with the following 1 3 steps: 

1 Calculate the constant elastic deflection vector {£/»}, using equation (7-3) with given 
{Fbo} and [CO 32 ], for all candidate contact points or gaps. 

2. Obtain the full constraint transformation matrix [Tj\ for all gaps from subroutine 
CNSTIF. The elements of sparse matrix [7/], (i.e., 0, -1, +1), are obtained 
immediately from the definitions of the various possible bearing support 
configurations. Matrix [T s ], used primarily in equations (7-18) and (7-19), is a subset 
of [7>] corresponding only to the trial set of assumed gap closures. 

3. Calculate the full vector of gap separation distances {Ag f } from equations (7-10) and 
(7-15) using matrix [7>], {Agf} relates to all gaps, opened and closed, whereas {Ag b } 
is the subset of {Ag f } corresponding only to the trial set of assumed gap closures. 

4. Assume the first trial set of assumed gap closures using subroutine CNCLSD. All 
radial gaps in the azimuthal plane nearest to the resultant applied radial force are 
initially assumed closed for both ball bearings and roller bearings. For the first pass 
with ball bearings, all axial gaps corresponding to the direction of the resultant axial 
load are assumed closed. Thus, for a positive axial load, gaps at positive outer ring 
faces are closed while gaps at negative faces remain open. 

5. Given the trial set of assumed gap closures, select the corresponding elements from 
matrix [7>] to form matrix [T,\ Using the same procedure, select the appropriate 
elements from vector {Ag f } to form vector {Ag b }. Thus matrix [T s ] and vector {Ag b } 
correspond to closed gaps only. 

6. Solve for the constraint forces {F cc } at each gap assumed closed using equation (7- 
19). Note that some of these constraint forces could be tensile (negative) if gaps are 
incorrectly assumed closed. This situation is corrected subsequently when a new set 
of assumed gap closures is defined. Note that equation (7-19) produces the exact 
linear solution for the given set of assumed gap closures. 

7. Expand the independent constraint force vector {F cc } into the full vector {F c } using 
equation (7-18). Note that {F cc } and {F c } could contain the effects of impossible 

tensile preload spring forces. 
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8. Calculate total gap separation distances {Ag} from equations (7-9) and (7-14) using 
matrix [T f ] with the matrix notation of equation (7-20). This {Ag} could contain the 
effects of impossible tensile preload spring forces. 

9. Calculate {Ag p } as the relative deflections at each preload gap obtained from the 
preliminary {Ag} calculation. Check for positive values of {Ag p } at preload gaps. If 
AgJk) > 0, calculate additional constraint forces {F cp } required to null the impossible 
tensile forces in the preload spring. These operations are performed in subroutine 
LFTOFF using equation (7-21). 

10. Recalculate {Ag} with the additional constraint forces {F cp } which null any preload 
spring tensile forces. This calculation uses the equation (8-1) below, rather than 

equation (7-20). 

{Ag} - m{U 3b } + [COssUFc + F cp }) + {G} (8-1) 

1 1 . Determine an updated trial set of gap closures based on two criteria: 

(a) A gap is closed if Ag{k) is less than or equal to a specified tolerance value 
which accommodates round-off error in calculating {Ag}. From equation (7- 
19), the constraint forces {F cc } are calculated to enforce Ag(k) = 0 at the gaps 
assumed closed. Because of round-off error, Ag(k) is usually calculated as a 
very small positive or negative number. The tolerance value is therefore a 
small positive number which may be machine dependent. 

(b) The trial set of gap closures determined from test (a) is modified to open any 
gap which has an independent constraint force F cc (k) less than zero. As noted 
previously, impossible tensile constraint forces can be generated by equation 
(7-19) if gaps are incorrectly assumed closed when they are actually open. 

These two tests attempt to define as closed those gaps which are in contact and 

which have contact compressive forces of zero or greater. 

12. Compare the updated set of assumed gap closures with the prior set of gap closures^ 
If the two sets of gap closures differ, return to step 5 and repeat the calculations wit 
the updated set of gap closures. When all F cc (k) > 0 and all Ag(k) > 0, continue to step 

13. 

13. When the closed gaps have been properly identified, calculate elastic deflections {U,b} 
using {Fbo} and, where applicable, (F e ) and {F cp } using equation (8-2). 

{ U,„} = [CO, 2 ] {Fto} + [CO, 3 ] {F cp + F cp } (8-2) 

Return to QSTATE to recalculate rolling element forces consistent with the kinematic 

constraint relations which now include the updated outer ring deflection vector {U,b}- 
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8.3 Iterative Solution Methodology 

The outermost iteration loop, as shown in Figure 8-2, determines overall elastic deflections of _the 
outer and inner ring configurations in response to applied loads or specified displacements. The 
iteration is managed by QSTATE, in conjunction with subroutine ERRCAL, which checks 
assumed elastic deflections against user specified limits and revises the displacements using a 
nonlinear solution methodology. 

The methodology incorporated into ERRCAL for the solution of the system of nonlinear 
equations uses a scheme from [9], This procedure is an acceleration method which does not use 
partial derivatives as with a Newton-Raphson approach. The calculation of partials was explicitly 
avoided due to the potential for large Jacobian matrices to be formed. For instance if a 
configuration with an outer ring of IBSCOR = 6 and 15 elements was used, the Jacobian would be 
22 x 15, containing 330 evaluations per iteration step. The procedure does calculate implicit 
partials, and employs a secant type estimate for the next root. 

This acceleration procedure attempts to minimize the difference between assumed and calculated 
elastic deflections at iteration n and n+J based on the following relationship. 

{*" +/ } - {**{*"})> (8 ' 3) 

which in the program, is implemented with the following algorithm 

= {6T>{^({x”»} +({/> - {ffms} ( g - 4 > 

where {/} is a unity vector and {#,} is a weighting factor, which for a single term is expressed as 

= l ( 8 - 5 ) 

' 1 - &,«*"}) 

where g„ = dgjdx t, is the partial derivative for the term under evaluation. In FEREBA, the actual 
value of any 9, is limited to 0.25, which was set to prohibit large changes in assumed values It 
was found that excessive corrections in iterates resulted in non-convergence for problems with 
many elastic displacement DOFs. 

During development testing of FEREBA, significant effort was placed in the selection of this 
elastic displacement iteration routine. Initially, it was theorized that setting the first loop values to 
zero, then calculating terms with DEFLOR and using these finite values for the second iteration 
would be sufficient. For linear ring configurations, this is valid, however, for the nonlinear 
IBSCOR types, the first set of displacements were “close” but not within a reasonable error 
amount. As such, the entire iteration problem can be stated as a refinement of the initial values 
determined in DEFLOR, and thus the limit on 6, to avoid large changes. 
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8.4 Variable List 
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FEREBA uses a storage scheme similar to SHABERTH, in that all of the primary problem 
variables are contained in several arrays. This section lists the contents of the arrays AVAR, 
BVAR, and BDAT, which in conjunction with the compliance matrices, store all of the significant 
program variables. Terms in these matrices can be printed for diagnostic purposes using the 
routine SYSCON, which is supplied with the source code. The first matrix listed is AVAR, which 
stores analysis variables. 


Variable 

AVAR ( 1 ) 
AVAR (2) 


AVAR (3) 

AVAR ( 4 ) 

AVAR (5-6) 
AVAR (7) 

AVAR (8) 

AVAR (9) 

AVAR (10) 
AVAR (11) 
AVAR (12) 
AVAR (13) 
AVAR (14) 
AVAR (15) 
AVAR (16) 
AVAR (17) 
AVAR (18) 
AVAR (19) 


Description 

Type of analysis requested, read in READAF 
NSLOT, type of quasi-static iteration run, set in READAF 
m 0 -> no iteration (normally not set) 

= 1 -> iteration without friction 
= 2 -> iteration with friction 

Number of requested elastic deflection iteration loops, read 
in READAF, used in QSTATE 

Required accuracy of elastic deflection iteration, read in 
READAF, used in QSTATE 

Number of requested inner ring force iteration loops, read 
in READAF, used in SLVBEQ 

Required accuracy of inner ring force iteration, read in 
READAF, used in SLVBEQ 

Flag to indicate if displacements (= 0) or loads (= 1) are 

specified, set in READAF *nair 

NPRINT, print flag for main analysis program, read in R 

IBSCOR, outer ring configuration (0, 9), read in READAF 

IBSCIR, inner ring configuration (0, 3), read in READAF 

MBO, outer ring force vector FBO size (defined in CODIM) 

MBI, inner ring force vector FBI size (defined in CIDIM) 

MC, outer ring constraint vector FC size (defined in CODIM) 
MCG, total number of outer ring gaps (defined in CODIM) 
Number of elastic deflection iteration loops, set in QSTATE 
Obtained accuracy in outer ring elastic deflection iteration 
Obtained accuracy in inner ring elastic deflection iteration 


The following 

AVAR (20) 
AVAR (21) 
AVAR (22) 


AVAR (23) 

AVAR (24) 

AVAR (25) 
AVAR (26) 
AVAR (27) 
AVAR (28-30) 


positions are control integers for rolling element quasi-static 
iteration (first five formerly stored in IBD in SHABERTH) 
IBRG, current bearing number being analyzed, set in BRGFOR 
ISTEP, current solution index , set in BRGFOR 

NARE, number of active rolling elements, set in BRGFOR & 
PREPAR (normally equal to total number of rolling elements 
except if cage only solution desired) 

IMAX, index of most heavily loaded rolling element with 
frictional solution, set in BRGFOR 

IRE, loop variable, current rolling element being analyzed 
(depends on ISTEP) , set in BRGFOR 
Number of iterations made in SLVBEQ 

Accuracy in force iteration obtained in SLVBEQ (ERRMAX) 
Number of times BEAREQ called by DNEQNF 
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The next matrix listed is BVAR, which stores specific rolling element variables. BVAR is a 3- 
dimension array. The first dimension is the pertinent term, as described below. The second 
dimension is the particular rolling element, maximum of 20. The third dimension is bearing 
number, maximum of 2. 


Element 


Description 


BVAR ( 1 ) 

BVAR (2) 

BVAR (3) 

BVAR (4) 

BVAR (5) 

BVAR (6) 

BVAR (7-8) 
BVAR (9-10) 

BVAR (11-12) 

BVAR (13-14) 


BVAR (15) 
BVAR (16) 
BVAR (17-20) 


BVAR (21-24) 


BVAR (25) 
BVAR (2 6) 
BVAR (27-30) 

BVAR (31-32) 

BVAR (33-34) 
BVAR (35-36) 


X-component of rolling element center position 
Y-component of rolling element center position 
X-component (axial) of rolling element rotational velocity 
Y-component of rolling element rotational velocity 
Z-component of rolling element rotational velocity 
Orbital speed of rolling element 
Outer/inner race-rolling element contact load 

For ball bearing, sine of outer/inner race rolling element 
contact angle. Zero for roller bearing. 

For ball bearing, cosine of outer/inner race rolling element 
, contact angle. Unity for roller bearing. 

For ball bearing, constant component of raceway groove 
curvature center distances ACON(l) and AC0N(2), calculated 
in BALLIN. Blank for roller bearing. 

UXP in center distance calculation for duplex bearing only 


Elastic deflections of master nodes, either axial or radial. 


Element 

Ball Bearing 

Roller Bearing 

17 

ORCC axial - UXORCC 

Left outer race radial - UROL 

18 

ORCC radial - URORCC 

Right outer race radial - UROR 

19 

IRCC axial - UXIRCC 

Left inner race radial - URIL 

20 

IRCC radial - URIRCC 

Right inner race radial - URIR 


calculated in DEFLOR & DEFLIR, checked in ERRCAL 


Outer/Inner ring rolling element forces, (axial or radial), 
depending on bear ing type (i refers to rolling elementjj 


Element 

*7 — 

Symbol 

2 J r : 

Ball Bearing 

Roller Bearing 

21 

FBO(i) 

Outer ring axial 

Left outer race radial 

22 

FBO (i+1 ) 

Outer ring radial 

Right outer race radial 

23 

FBI (i) 

Inner ring axial 

Left inner race radial 

24 

FBO (i+1) 

Inner ring radial 

Right inner race radial 

. / t— « t'i -r _• _ ptt mrrs 


calculated in CALFOR/ SUMFIR, stored in FBO/FBI in SLVBEQ . 

FW (gyro force) for outer ring, (not used) , refed in SUMFIR 
ZM for outer ring, calculated in ROLLEQ, used in CALFOR 


Friction force components acting on inner ring; FW, FZ, YM, 
ZM, used in SUMFIR (refed but not used) 

EHD film thickness at most heavily loaded outer/inner 
raceway/slice, calculated in FMIX/FMIXR (not used) 
Outer/inner raceway contact area, used in FMIX and FMIXR 
PMAX - in BALLEQ, maximum hertz stress, and in FMIXR, max 
hertz stress at most heavily loaded slice 



The main storage of problem variables is in matrix BDAT, as listed below 


Variable 

BDAT ( 1 ) 
BDAT (2) 
BDAT (3) 
BDAT ( 4 ) 
BDAT (5) 
BDAT (6) 
BDAT (7) 
BDAT ( 8 ) 


Description 

Bearing type (0, 1, 2, for ball, roller, or duplex) 
Orientation angle of first rolling element (always 0) 
Rotation speed of inner ring/shaft 
Rolling element diameter 

Pitch diameter to rolling element center 
Number of rolling elements 
Diametral clearance 

Initial unloaded contact angle for ball bearing, based on 
diametral clearance, calculated in CALCON 


BDAT (9) 

BDAT (10) Roller length end to end 

BDAT (11) Roller end sphere radius 

BDAT (12) Roller included angle 

. The next six items are not used in quasi-static 
retained for use in future frictional analyses 


code, but the positions are 


BDAT (13) Outer race CLA surface roughness 

BDAT (14) Inner race CLA surface roughness 

BDAT (15) Rolling element CLA surface roughness 

BDAT (16) Outer race asperity slope 

BDAT (17) Inner race asperity slope 

BDAT (18) Rolling element asperity slope 


BDAT (19) 
BDAT (20) 
BDAT (21) 
BDAT (22) 

BDAT (23) 
BDAT (24) 
BDAT (25) 
BDAT (26) 
BDAT (27-30) 
BDAT (31) 
BDAT (32) 
BDAT (33) 
BDAT (34) 
BDAT (35) 

The next six 
retained for 
BDAT (36) 

BDAT (37) 
BDAT (38) 
BDAT (39) 
BDAT (40) 
BDAT (41) 


Roller crown radius 

Roller flat length or ball bearing outer race curvature 

Ball bearing inner race curvature or number o 

roller/raceway axial slices (maximum of 20, set in READAF) 

Roller outer raceway effective length 

Roller inner raceway effective length 

Roller outer raceway crown radius 

Roller inner raceway crown radius 

Cylindrical roller bearing outer ring flange angle 

Cylindrical roller bearing inner ring flange angle 

Cylindrical roller bearing outer ring end play 

Cylindrical roller bearing inner ring end play 

Cylindrical roller bearing flange inversion flag 
items are not used in quasi-static code, but the positions are 

use in future frictional analyses 

Cage type (-1, 0, +1) for outer ring land, inner ring land, 

or rolling element riding 

Cage rail land diameter 

Cage single rail width 

Cage rail land diametral clearance 

Rolling element cage pocket diametral clearance 

Cage weight 
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Contents of main problem variable storage matrix BDAT, continued 


Variable 

BDAT (42-43) 

BDAT (44-45) 
BDAT (4 6-50) 

BDAT (51) 
BDAT ( 52 ) 
BDAT (53) 
BDAT (54) 
BDAT (55) 
BDAT (56) 
BDAT (57) 
BDAT (58-65) 
BDAT (66-70) 

BDAT (71-75) 

BDAT (76-80) 

BDAT (81-85) 


Description 

Input (42) and output (43) values of cage degrees of 
freedom, set in PREPAR 

Life multipliers for outer/inner ring 

Input values of applied displacements (x-trans, y-trans, z 
trans, y-rot, z-rot) 

Inner ring mean outer diameter 

Outer ring mean inner diameter 

Outer ring mean outer diameter 

Carrier mean inner diameter 

Axial gap between outer ring and carrier 

Initial axial distance between duplex ball bearing centers 
Initial axial gap between duplex bearing outer rings 

Elastic Moduli for shaft, inner ring, rolling elements, 
outer ring, and housing 

Poisson's ratio for shaft, inner ring, rolling elements, 
outer ring, and housing 

Density for shaft, inner ring, rolling elements, outer ring, 
and housing 

Thermal Expansion Coefficient for shaft, inner ring, rolling 
elements, outer ring, and housing 


BDAT (86-87) ^ M M . 

BDAT (88-92) Input values of applied forces (Fx, Fy, Fz, My, mz) 

The next • eighteen items are not used in quasi-static code, but the positions 
are retained for use in future frictional analyses 

BDAT (93) Lubricant code - if = 0, implies no friction 

= +1 i-o ±9 lubricant codes (friction) 


BDAT (94) 

BDAT (95) 

BDAT (96) 

BDAT ( 97 ) 

BDAT (98) 

BDAT (99-101) 
BDAT (102) 
BDAT (103) 
BDAT (104) 
BDAT (105) 
BDAT (106) 
BDAT (107) 
BDAT (108) 
BDAT (109) 
BDAT (110) 


= ±10 implies dry friction 
the ± refers to traction mod (see SHABERTH manual) 
Kinematic viscosity at 100° F 
Kinematic viscosity at 210° F 
Density at 15.5° C 

Lubricant coefficient of thermal expansion 
Lubricant thermal conductivity 
Calculated data from LUPROP and LUBCON 

Outer raceway lubricant replenishment layer thickness 
Inner raceway lubricant replenishment layer thickness 
Percent lubricant in bearing cavity 
BDAT (104) /100, also XCAV in BCON & CRCON 
Raceway-Rolling element asperity friction coefficient 
Cage friction coefficient 

Outer race flange replenishment layer thickness 
Inner race flange replenishment layer thickness 
Flange-Roller end friction coefficient 


BDAT (111) 

BDAT (112-113) 
BDAT (114) 

BDAT (115) 


Inner ring speed in rad/sec, calculated in CALCON 

Rolling element y-axis inertia, calculated in CRCON 
Rolling element x-axis inertia, calculated in BCON 
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Contents of main problem variable storage matrix BDAT, continued 

Variable Description 

BDAT (116) Rolling element mass, calculated in BCON 

BDAT (117) Sine of initial contact angle BDAT(8) calculated in CALCON 

BDAT (118) Total curvature of ball bearing (B) , calculated in CALCON, 

used for contact angle calculation and in PRELOD 

BDAT (119) 

BDAT (120-121) Storage for variable EMOD in SHABERTH s/r CONS, calculated 

in CALCON, used in FMIX et al routines 
BDAT (122) RMS values of surface roughness, array RMS in CRCON 

BDAT (123) Sum of BDAT (127) and BDAT (128), variable CRVS in BCON (= A 

or BD) , RMS (2) in CRCON 

BDAT (124) Cosine of initial contact angle BDAT(8) calculated in CALCON 

BDAT (125-126) Sine of roller bearing flange angles, array SINF in CRCON 

BDAT (127-128) Osculation parameters, array CRV in BCON; Cosine of roller 

bearing flange angles, array COSF in CRCON 
BDAT (129-130) Array ZEMA in BCON (See Jones 13.189 & 13.190); Flange 

height parameter, array ARMS in CRCON 

BDAT ( 131-132 ) Raceway/rolling element contact moduli calculated in CALCON 
BDAT (133) . Variable Z10 defined in CALCON 

BDAT (134) Variable GP = BDAT (4) / BDAT (5) calculated in BCON & CRCON 

BDAT (135-136) 

BDAT (137) Variable Y calculated in BCON (initial IRCC radius, eqn 

13.177 in Jones); WN2 = BDAT (22) /2 in READAF 
BDAT (138) Ball drag constant Zl, used in BALLEQ 

BDAT (139) Ball drag constant Z2, calculated in BCON & CRCON 

BDAT (140) Ball drag constant Z3, used in BALLEQ 

BDAT (141-142) Array W in CRCON, slice width 

BDAT (143-144 ) Ratio involving raceway curvatures, array RDEF in BCON 
BDAT (145-146) Ratio involving raceway curvatures, array REX in BCON; 

raceway array C2 in CRCON 

Delta angular position of each rolling element 
152) Contact deformation values CDEF(3,2), calculated in ABDEL 

155) Array RX in CRCON 

Variable PREF in BCON, equal to .05 * rolling element dia 
Relative inner ring speed (WREF) , calculated in CALCON 
159) Variable TF & TM in CALCON (see BRGAX for more details) 
Variable TM in CALCON 

Variable WREFYZ in BCON & CRCON, .2 * WREF (BDAT (157)) 

163) Equated to BSC(l-3) in SPRING 

Variable RPZ in CALCON, used in GUESCG 

167) 

Variable CRP in GUESCG, could be number of cage DOF 
Variable DELYZ in GUESCG, vector sum of cage radial 
displacements 

BDAT (170) Variable THETA in GUESCG, angle of cage radial displacement 

BDAT ( 171-175) Calculated 5-axis inner ring displacements (x-trans, y- 

trans, z-trans, y-rot, z-rot) 

BDAT (17 6-180) Calculated 5-axis inner ring forces (Fx, Fy, Fz, My, Mz) 


BDAT (147) 
BDAT (14 8 — 
BDAT (153) 
BDAT (154- 
BDAT (156) 
BDAT (157) 
BDAT (158— 
BDAT (159) 
BDAT (160) 
BDAT (161- 
BDAT (164) 
BDAT (165) 
BDAT (166— 
BDAT (168) 
BDAT (169) 
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Contents of main problem variable storage matrix BDAT, continued 


Variable 

BDAT (181-184) 
BDAT (185) 

BDAT (186) 

BDAT (187) 


BDAT (188-189) 
BDAT (190) 


Description 


Change in bearing diametral clearance due to temperature, 
calculated in TMPFIT, used in CALCON, GROLLB, ROLLIN 
Operating diametral clearance = BDAT (7) + BDAT (185) [not 

used] 

Change in clearance between outer ring OD and housing for 
IBSCOR configurations with deadband, calculated in TMPFIT, 
used in GAPDEF 

force, used only with duplex 


axial preload 

preload spring stiffness. 


used with ball 


Specified 
bearings 

BDAT (191) Specified axial 

bearings 

■BDAT (192) Resulting elastic axial displacements due to preload (UXP ), 

used only with duplex bearings 

BDAT (193-200) 

The next thirty items are not used in quasi-static code, but the positions are 

retained for use in future frictional analyses 

Maximum EHD film thickness HMAX defined in FMIXR 
Miniscus distance variable DIST defined in FMIX & FMIXR 
EHD reduction starvation factor PHIS defined in FMIX & FMIXR 
EHD reduction thermal factor PHIT defined in FMIX & FMIXR 
Outer/inner raceway lubrication life factor 
Outer/inner raceway and overall bearing fatigue life 
Equated to BDAT (215) in GUESCG 

Variable WCAGE in GUESCG, angular speed of cage 


BDAT (201-202) 
BDAT (203-204) 
BDAT (205-206) 
BDAT (207-208) 
BDAT (209-210) 
BDAT (211-213) 
BDAT (214) 

BDAT (215) 

BDAT (216-218) 
BDAT (219-220) 
BDAT (221-222) 
BDAT (223-224) 
BDAT (225) 

BDAT (226) 

BDAT (227) 

BDAT (228) 

BDAT (229) 

BDAT (230) 


Array XCGF in CRCON 

Composite raceway surface roughness, calculated in CALCON 
Composite raceway asperity slope, calculated in CALCON 
Outer ring flange CLA surface roughness (array SIGF in 
BCON/CRCON) 

Inner ring flange CLA surface roughness 
Rolling element end CLA surface roughness 

Outer ring flange asperity slope (array SGF in BCON/CRCON) 
Inner ring flange asperity slope 
Rolling element end asperity slope 


BDAT (231-232) 
BDAT (233-234) 
BDAT (235-236) 
BDAT (238) 

BDAT (239-240) 
BDAT (241-242) 
BDAT (243-244) 
BDAT (245-246) 
BDAT (247-248) 
BDAT (249-270) 


Array HT in CRCON 
Array RACE in CRCON 
Array RCFY in CRCON 

Variable XLCG in CRCON (half roller length) 

Array RRFY in CRCON 
Array RRF in CRCON 
Array RCFX in CRCON 

Array GAP in SETGAP, roller end gap 
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Contents of main problem variable storage matrix BDAT, continued 


Variable 

The next thirty 
retained for use 
BDAT (271-274) 

BDAT (275-278) 

BDAT (279-282) 

BDAT (283-286) 

BDAT (287-294) 
BDAT (295-298) 
BDAT (299-300) 


BDAT (301-330) 
BDAT (331-360) 


Description 


items are not used in quasi-static code, but the positions are 
in future frictional analyses 

Variable FLFAC(l) in ROLLEQ passed to FLMIX as DISTH, 
miniscus distance variable 

Variable FLFAC (5) in ROLLEQ passed to FLMIX as PHISH, EHD 
film reduction starvation factor 

Variable FLFAC (9) in ROLLEQ passed to FLMIX as PHITH, EHD 
film reduction thermal facor 

Variable FLFAC (13) in ROLLEQ passed to FLMIX as RATIOH, 
QASP/QTOT for roller bearing 


Flange/roller conductivity terms 
Outer/inner film thickness to 
calculated in LIFE 


surface roughness ratio. 


Sine of rolling element angular position, set in CALCON 
Cosine of rolling element angular position, set in CALCON 


In the next six items, the first 20 positions refer to 
outer race, the last 20 positions refer to the inner race 

_ . . ~ . . i i _ 1 J _ _ — — up 


the roller bearing 


BDAT (361-400) 
BDAT (401-440) 
BDAT (441-480) 
BDAT (481-520) 
BDAT (521-560) 
BDAT (561-600) 


Array DK in SLICES, roller slice radius 
Roller bearing crown drops CD, calculated in SLICES 
Slice half width constant BMI, calculated in SLICES 
Array HM in SLICES, called from CRCON 

Product of terms used to calculate film thickness (not used) 
Array RK in SLICES 


BDAT (601-605) 


BDAT (606) 

BDAT (607-650) 
BDAT (651-690) 

BDAT (691-750) 


Temperatures of shaft, inner ring, rolling elements, outer 
ring, and carrier/housing. Note that flange and lube 
temperatures used in SHABERTH are not defined here. 
SHABERTH initially stores temperatures in an array TB, then 
transfers to these locations in routine SHABE. 

Reference temperature for fit calculations, set to 70°F in 
CALCON, used in TMPFIT 

Variable HZ in ROLLEQ, hertz contact stress, first 20 are 
for outer raceway, second 20 are for inner raceway 
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ENCLOSURE (1) 


STATIC LOAD -DISPLACEMENT TEST OP SINGLE BALL BEARING 
D. H. Merchant 23 July 1994 


TEST DESCRIPTION AND RESULTS : 

A discussion of the special relations between forces 
applied to the test rig and the statically equivalent forces 
and moments input to FEREBA is provided in Attachment 1. 
Attachment 2 is a listing of the FEREBA input data file used 
for the pretest analysis for the no-deadband configuration. 
The axial preload spring stiffness of 73,400 lb/in 
represents the slope of the nonlinear load-displacement 
curve for small loads. This Belleville spring bottoms with 
a maximum load of about 3800 lb at a maximum stroke ("Outer 
Ring/Hsg Axial Gap") of 0.093 in. For the deadband test 
configuration, the carrier inner diameter exceeds the outer 
ring outer diameter by 0.0002 in; this results in a measured 
radial deadband of 0.0001 in. 

Attachment 3 comprises the informal notes, for both 
bearing tests, prepared by the Test Engineer. Raw data 
consisting of applied forces, measured by four load cells, 
and corresponding displacements, measured by four Bently 
probes, are listed on page 3-4 for Test 1 and on page 3-7 
for Test 2. The XF-YF-ZF coordinates for the transducer 
locations listed in Table 1 apply to the FEREBA analysis. 
To order the test rig locations for consistency with those 
for the FEREBA analysis, it is necessary to interchange the 
values for B2 and B3 and for L2 and L3. 

As noted in Attachment 1, the forces applied to the 
test rig are the carrier support reactions in the FEREBA 
analysis. The FEREBA input loads must thereforce be 
calculated from the applied test rig loads (plus distributed 
carrier weight) by appropriate transformation relations. 
These equations are coded, for example, in Attachment 4a for 
Case 6 of the deadband configuration (Test 2) . The 
corresponding forces and moments input to FEREBA for this 
case are listed in Attachment 4b. A complete set of test- 
rig applied forces and statically equivalent FEREBA forces 
and moments are listed in Tables 2 and 3 for both tests. 
Forces for the five runs for each load case are averaged. 
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For comparison with FEREBA calculated displacements, 
the test rig displacements must be converted to equivalent 
translations and rotations relative to the shaft inner ring. 
This conversion is accomplished by the transformation 
equations coded, for example, in Attachment 5a for Case 1 of 
the no-deadband configuration (Test 1) . The corresponding 
equivalent shaft displacements for the case are listed in 
Attachment 5b. A complete set of these test rig "measured" 
displacements are listed in Tables 4 and 5 for both tests. 
Displacements for the five runs for each load case are 
averaged . 

The displacements measured on the test rig by the 
Bently probes do not include the initial axial displacement 
due to the carrier weight of 24 lb. For comparison with 
test rig "measured" displacements, the axial displacement 
to preload (Ux = 0.0006 in) must be subtracted from the 
FEREBA displacements calculated with applied loads plus 24- 
lb axial preload. The resulting relative shaft displacements 
are listed in Tables 4 and 5 as FEREBA output displacements. 


CONCLUSIONS : 

1. The forces applied to the test rig for both tests 
show very small run- to- run variations from the intended 
forces . 

2. The displacements measured for the no-deadband 
configuration show acceptably small run- to- run variations. 
In particular, the four runs without applied load (Case 0) 
show excellent repeatability. 

3. The displacements measured for the 0.0001 -inch 
deadband configuration show rather large run- to- run 
variations. In particular, the six runs without applied 
load (Case 0) show poor repeatability. This probably 
indicates stiction in the bearing load path. 

4 . One variable which could not be accounted for 
between test and analysis is the alignment of radial force 
relative to ball locations. The potential magnitude of this 
effect was not evaluated. 

5 . FEREBA did not converge to input forces and moments 
for several cases. Correlation between test "measured" 
displacements and FEREBA output displacements is difficult 
to assess because of the relative errors in translation and 
rotational displacements. 
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6. An alternate procedure is available with FEREBA to 
correlate measured and calculated data. It is possible to 
input "measured" displacements to FEREBA and calculate 
corresponding shaft forces to compare with test values 
designated as FEREBA input forces. This procedure would 
avoid the previously noted convergence problems since no 
convergence is involved with the displacement input option. 
Whether assessing correlation becomes easier with this 
procedure was not determined. 
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ATTACHMENT 1 


REPRESENTATION OP TEST RIG FOR FEREBA PRETEST ANALYSIS 
Version 1 D. H. Merchant 12-15-93 


COORDINATE SYSTEM A ND CONSTRAINT^ : 

In FEREBA, the thrust force is applied in the + X 
direction to the inner ring and shaft at the center of the 
unmounted bearing. This corresponds to the + Z d ^ection in 
the ANSYS finite -element model. FEREBA then automatically 
transforms the ANSYS reduced stiffness matrix ^^y^ndrical 
coordinates from the ANSYS coordinate system J^® A -ZA) to 
the required FEREBA cylindrical coordinate system (XF-RF- 
*F? The constraints applied to the ANSYS carrier model 
must be consistent with the FEREBA coordinate system. 

The solution procedure in FEREBA is predicated on a 
bearing system which has the carrier fixed m inertial space 
and which has forces or displacements prescribed at the 
inner ring and shaft. Whether forces or displacements are 
prescribed, the inner ring/shaft undergoes five independent 
snail displacements to result in rolling element forces^ 
For the test rig, the inner ring and sha : £t * re ™ 

inertial space and forces are applied to the carrier to 
result in rolling element forces. 

To represent the bearing and carrier of the test rig 
using the FEREBA solution procedure, it is necessary to 

replace the forces applied to the test ANSYS /FEREBA 

reaction forces at the carrier supports of the ^ sy S/FEREBA 
model. By constraining six degrees of freedom on the ANSYS 
carrier in a statically determinate manner and by applying 
lopropriate forces and moments at the center of the 

unmounted bearing, the reaction f< orce . s ' £ e * ] ^h^tlstTrig 

can be made identical to the forces applied to the test rig 
carrier. Four of the six constrained degrees of freedom in 
the ANSYS model must be exactly those at which the three 
axial forces and one radial force are applied to the 
rig carrier. And the forces and moments applied at the 
center of the FEREBA bearing must exactly equilibrate the 
forces applied to the test rig carrier. 

The deflections measured by the test rig are relative 
deflections between the fixed inner ring and the movable 
carrier. The measured deflections are therefore equivalent 
to those calculated by FEREBA so long as there is negligible 
elastic 6 deflection of the test rig's shaft connecting the 
inner rinq to ground. Since this shaft is solid steel with 
diameter Exceeding that of the inner ring, the shaft elastic 
deflections are justifiably neglected. 
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ATTACHMENT 2 


FEREBA INPUT DATA FOR PRETEST ANALYSIS (CASE 6) 


It f 9 fell fearing, IlSCOft/IISCIt ■ 3/1, Fa ■ 828, Fr • 600 


Oatt this analysis fila written: 
TIae this analysis file written: 
f ester, asc 


Thu ft b 3, 199* 
10:36:50 

ferabfr.ssc 


CONTROL 

1 

0 

0 

50.00000 
0.001000000 

100.0000 

1 .0000006*04 
2.000000 

3.000000 
0.00000006*00 

GENERAL 

0.00000006*00 

3.366000 

14. 00000 

1.000000 
1.000000 
2.974000 
3.784000 
4.331000 
4.331000 

9.30000026*02 

73400.00 

SAIL MG 

0.6250000 
2.00000016*03 
0.5300000 
0.5200000 
NATL PROP 
440C 
52100 
52100 
52100 
52100 

3.00000006*07 
3.00000006*07 
3.00000006*07 
3.00000006*07 
3.00000006*07 
0.3000000 
0.3000000 
0.3000000 
0.3000000 
0.3000000 
0.2830000 
0.2830000 
0.2830000 
0.2830000 
0.2830000 
8.00000006*06 
8.00000006*06 
8.00000006*06 
8.00000006-06 
8.00000006*06 
TEMP DATA 

70.00000 
70.00000 
70.00000 
70.00000 

70.00000 

LOADS 

828.0000 
0.00000006*00 
600.00000 
1500.000 
0.00000006*00 


Analysis Type 

Miring Typt (0, 1, or 2) 

Analysis Units 

Ring Dfsp Iteration loops 

ling OUp Iteration Error 00 

force Iteration Loops 

force Iteration Error PQ 

Analysis Print flag 

Outer ling Configuration 

Inner ting Configuration 

Shaft Speed 

Maring Pitch Diaaeter 
Suiter of Rolling Eleaenta 
Outer ting Life Multiplier 
Inner ting Life Multiplier 
Inner ting Outer Ofaaettr 
Outer ting Inner Ofaaettr 
Outer ting Outer Ofaaettr 
Cerrier/Heg Inner Ofaaettr 
Outer ting/Hsg Axial Gap 
Preload Spring Stiffness 

Mil Diameter 
Oiaaetral Clearance 
Outer Raceway Curvature 
Inner Raceway Curvature 

Shaft Material 
Inner ting Material 
Elaaents Material 
Outer ting Material 
Car/Nang Matarial 
Shaft Elastic Modulus 
Inner ting Elastic' Modulus 
Eleaenta Elastic Modulus 
Outer ting Elastic Modulus 
Car/Nsng Elastic Modulus 
Shaft Poisson Ratio 
Inner ting Poisson Ratio 
Elsaents Poisson Ratio 
Outer ting Poisson Ratio 
Car /Rang Poisson Ratio 
Shaft Density 
Inner Ring Density 
Eleaents Density 
Outer Ring Density 
Car/Nsng Density 
Shsft Thenssl Exp Coeff 
Inner Ring The real Exp Coeff 
Elsaents Thersal Exp Coeff 
Outer Ring Thersal Exp Coeff 
Car/Nsng Thermal Exp Coeff 

Shaft Temperature 
Inner Ring Temperature 
Eleaents Temperature 
Outer Ring Tenpersturt 
Car/Nsng Tespersturs 

X*Axis (Thrust) force 
Y-Axis force 
Z-Axis force 

XZ-plane Moment (About Y) 

XT -plane Moment (About 2) 
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OaoiNAL PAGE fS 

or poor quality 




4. JL A 


BEARING TEST NOTES - ~GK WONo 

The two bearing teste were conducted .t A rone in Ed Mack's lab. The Tirst 
one on 6-17-94, the second on 6-21-94. 

Calibration of load cells by using -“age -bstition 
with exitation adjusted to lO.OOvdc. 

Calibration of Bently probes by d PPlV ‘"9 various thickness ot 
record output voltages. Shims used : 0.010 , 0.015 , 0.0-0 » ’ 

0.048" and 0.059". 

Test data were taien 5 times on each load condition (case) . Loads were 
1800 lbf due to cable yielding at about aAO 1 

• M hant and fin t - - y ^reviewed the raw Bently voltage outputs on each 
ctse aUer Jhe itt bearing test. It was determined that the data seemed to 
consistent enough that we will go into the 2nd bearing test. 

. - tz \ /-iilt LN° was used to freeze the center 

h sez v 

The Tnd bearing was installed after freezing the mounting shaft in LN2 for 
S minutes. The bearing elided Into place with no apparent effort. 

interference in .'fi* thicker’ shims. This Approach seemed to improve reading 

T'Frfo -: UVn ^“g S -1 SilriZlrr^ti^yW^S^* 

. i'Ji: .Si^^-tT^STnStSKS ?d th the 

| 1st bearing cal, then process the test data. 

’ There ‘U‘ 

There for the thin target discs were used. 


The following is 


the brief description of files generated: 


READ. ME - This file. 

BENTCAL.WQ1 - The Quart tro file of 


th< 


2 bearing test's Bently probe cals. 


RAWBEAR1 . DAT 
BEAR I NG 1 . DAT 

L I NBEAR 1 . PRN 


PC 

II 


for bearing test i. 


case sin 


- Raw data taken on . - ■— - . . 

- Raw data “clean up" to just one data entry per load 
data is very stable during each case load. Test 1. 

- Linearized data after ran through Russ Miller sBENT_LIN.bas 
program using calibration file BCAL.cal. Test 1. 

BEARING^WQl^-'^Quattro* f iTe^wi th LINBEAR1.PRN impor ted with headings. Test 
IcAL^cill - Bently cal. Test 2 with the last two thickness shim output 
BCAL-.C reading ..normalized" from the first calibration. 

OWGINAL PAGE fS 
OF POOR QUALITY 
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— — . . ^ — t-T^re - r . ... _ 

Patn: u: \Qv*+*-> 

Files BCAL1 .CAL 306 .a.. 


6-21-94 


Page 1 


5,7 

0.01 , 0.015 

5.61 , 6.70 

5 . 47 , 6 . 65 

4.45 , 5.56 

1.56 , 2.37 


0.02 , 
7.96 , 
7.90 , 
6.83 , 
3. 44 , 


0.029 , 
10.03 , 
10.26 , 
9.37 , 
5.50 , 


0 

• 

039 

1 

2 

. 12 

1 


. 26 

1 

i 

.42 


7 

.41 


0.048 , 
13.89 , 
1 4 . 00 , 
12.96 , 
9.26 , 


0.0 

15. 

15. 

14. 
11.38 


ORIGINAL PAGE fS 
OF POOR QUALITY 


cn ~j u cn 


Files LINBEAR1 . PRN 2,551 .a.. 6-22-94 


Page 1 


57613.4 , 

50178.7 , 

58692.9 , 

59503.0 , 

62755.1 , 

63128.2 , 

63377.7 , 

63954.9 , 

64485.8 , 

65201.4 , 

66303.3 , 

66790.8 , 
66917.1 , 

67328.5 , 

68044.3 , 

68425.8 , 

69455.8 , 

69764.4 , 

70216.0 , 
7<_>426 . O , 
70754 .3 , 

711 45 .6 , 

71434.0 , 

71572.7 , 

71837.0 , 

72361.7 , 

72795.8 , 
73355.2 , 

73831.5 , 

73961-0 » 


0.0293 , 
0.0241 , 

0.0213 , 
0.0457 , 
0.0230 , 
0.0197 , 
0.0292 , 
0.0236 , 
0.0207 , 
0.0459 , 
0.0229 , 

0.0198 , 
0.0291 , 

0.0233 , 
0.0204 , 
0.0452 , 
0.0231 , 
0.0197 , 
0.0237 , 
o . o 2 0 1 , 
0.0453 , 
0.0232 , 
0 .0196 , 
0.0292 , 
0.0240 , 
0.0212 , 
0.0456 , 
0.0238 , 
0.0201 , 
0.0292 , 


0.0261 , 
0.0287 , 
0.0296 , 

0 . 022 1 , 

0 . 0257 , 
0.0272 , 

0 . 0280 , 
0.0289 , 

0 . 0299 , 
0.0220 , 
0.0257 , 
0.0272 , 
0.0281 , 
0.0293 , 
0 . 0302 , 
0.0221 , 
0.0258 , 
0.0273 , 
0 . 029 i , 
0.0301 , 

0.0221 , 
0.0257 , 
0.0273 , 
0.0281 , 
0.0289 , 
0.0298 , 
0.0221 , 
0.0253 , 
0.0271 , 

0.0280 , 


0.0281 , 
0.0328 , 
0.0341 , 

0.0227 , 
0.0291 , 

0.0317 , 
0.0313 , 
0.0335 
0.0347 , 
0.0226 , 
0.0293 , 
0.0316 , 
0.0317 , 
0.0335 , 
0.0349 , 
0.0227 , 
0.0289 , 
0.0316 , 
0.0332 , 
0.0345 , 
0.0227 , 
0.0289 , 
0.0317 , 
0.0316 , 
0.0330 , 
0.0342 , 
0.0227 , 
0.0287 , 

0.0312 , 
0.0316 , 


0.0340 , 
0.0267 , 
0.0237 , 
0.0502 , 
0.0280 , 
0.0238 , 
0.0318 , 
0.0260 , 
0.0231 , 

0.0505 , 
0.0279 , 
0.0240 , 
0.0315 , 
0.0258 , 
0.0226 , 
0.0499 , 
0.0282 , 
0.0239 , 
0.0261 , 
0.0230 , 
0.0499 , 
0.0282 , 
0.0238 , 
0.0316 , 
0.0265 , 
0.0235 , 
0.0501 , 
0.0289 , 
0.0244 , 
0.0317 , 


62 , 

59 , 
61 , 

600 , 
301 , 

301 , 

0 , 
61 , 

60 , 

599 , 

30 1 , 

300 , 

0 , 
61 , 
61 , 

600 , 

300 , 

301 , 
60 , 

59 , 
602 , 

300 , 

30 1 , 
0 , 

61 , 

60 , 
599 , 

30 1 , 

302 , 
0 , 


61 , 

61 

» 


59 , 

59 


60 

61 , 

61 


18-: 

116 , 

116 

» 

6*: 

300 , 

302 

T 

6'". 

301 , 

301 

1 

ISC 

0 , 

0 

f 


60 , 

59 


6C 

59 , 

59 

1 

18C 

114 , 

114 

1 

6C 

301 , 

300 

1 

5- 

300 , 

301 

1 

ieo 

0 , 

0 

1 


60 , 

59 

1 

6C 

60 , 

60 


1 7 C 

113 , 

114 

1 

6C 

299 , 

299 

1 

6C 

302 , 

301 

J 

18C 

60 , 

59 

1 

60 

58 , 

60 

1 

iec 

1 15 , 

114 

f 

60 

300 , 

301 

1 

6‘. 

301 , 

301 

1 

1 SO 

0 , 

0 

1 


60 , 

60 

1 

60 

60 , 

61 

* 

160 

115 , 

1 15 

1 

60 

298 , 

298 

1 

5 C 

299 , 

301 

1 

180 

0 , 

0 

1 
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ORIGINAL PAGE ?S 
OF POOR QUALITY 


Alo DEADBAA/P 

BEARING TEST NO 1 @ A ZONE 6-1 7-94 
(5 runs were performed on each case ) 

B1 ...B4 converted to INCH 
L1....L4 in LBf 



\/iew 
p owa/ 






L — 1 > /c a</ c « 

i\s 







£ -p def I* «- 


J <+J 


TiME 

61 

B2 

B3 

B4 

LI 

L2 

L3 

L4 

CAGE 

53643-4— 0 0293 Cr0261 - 

-00291 

00340 62 - 



— 61 

6- 


58178.7 

0.0241 

0.0287 

0.0328 

0 0267 

59 

59 

59 

600 

5 

58692.9 

0.0213 

0.0296 

0 0341 

0 0237 

61 

61 

61 

1801 

4 

59503 0 

0.0457 

0.0221 

0.0227 

0 0502 

600 

116 

116 

602 

6 

62755.1 

0 0230 

0 0257 

0 0291 

0.0280 

301 

300 

302 

600 

2 

63128.2 

0.0197 

0.0272 

0 0317 

0.0238 

301 

301 

301 

1803 

1 

63377.7 

0 0292 

0.0280 

0.0313 

0.0318 

0 

0 

0 

0 

0 

63954.9 

0.0236 

0:0289 

0.0335 

0.0260 

61 

60 

59 

601 

5 

64485 8 

0.0207 

0.0299 

0.0347 

0.0231 

60 

59 

59 

1802 

4 

65201 .4 

0.0459 

0.0220 

0.0226 

0.0505 

599 

114 

114 

604 

6 

66303 3 

0.0229 

0.0257 

0.0293 

0.0279 

301 

301 

300 

599 

C. 

66790 8 

0.0198 

0 0272 

0 0316 

0 0240 

300 

300 

301 

1800 

1 

66917.1 

0.0291 

0.0281 

0 0317 

0.0315 

0 

0 

0 

0 

£• 

673285 

0.0233 

0.0293 

0.0335 

0 0258 

61 

60 

59 

601 

5 

68044.3 

0 0204 

0 0302 

0.0349 

0.0226 

61 

60 

60 

1798 

4 

68425.8 

0 0452 

00221 

0.0227 

0.0499 

600 

113 

114 

601 

6 

69455 8 

0.0231 

0.0258 

0.0289 

0.0282 

300 

299 

299 

600 

2 

69764 4 

0 0197 

0.0273 

0.0316 

0.0239 

301 

302 

301 

1802 

1 

70216 0 

0.0237 

0 0291 

0.0332 

0.0261 

60 

60 

59 

600 

5 

70426.0 

0.0207 

0 0301 

0.0345 

0.0230 

59 

58 

60 

1801 

4 

70754 3 

0 0453 

0 0221 

0.0227 

0.0499 

602 

115 

114 

600 

6 

71145.6 

0 0232 

0.0257 

0 0289 

0.0282 

300 

300 

301 

602 

2 

71434.0 

0.0196 

0.0273 

0.0317 

0.0238 

301 

301 

301 

1800 

1 

71572.7 

0.0292 

0 0281 

0.0316 

0.0316 

0 

0 

0 

0 

0 

71837.0 

0.0240 

0.0289 

0.0330 

0.0265 

61 

60 

60 

600 

5 

72361 .7 

0.0212 

0.0293 

0 0342 

0.0235 

60 

60 

61 

1801 

4 

72795 8 

0 0456 

0 0221 

0 0227 

0.0501 

599 

115 

115 

601 

6 

73355 2 

0.0233 

0 0253 

0.0287 

0.0289 

301 

298 

298 

599 

2 

73831 .5 

0 0201 

0.0271 

0.0312 

0.0244 

302 

299 

301 

1800 

1 

73961.0 

0 0292 

0.0280 

0 0316 

0.0317 

0 

0 

0 

0 

0 
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File; ECAL2 -CAL 


6 - 22-94 


Page 1 


.01 , 0.015 , 

.87 , 6.71 , 

.53 , 6.60 , 

.45 , 5.56 , 

1.69 , 2.61 , 


0.02 , 
7.85 , 
7.80 , 
7.09 , 
3.61 , 


0.029 , 

0.039 


10.02 , 

11.92 


9.96 , 

11.70 


9.08 , 

10 . 77 


5 . 18 , 

6.98 



0.048 

1 

0.059 

13.69 

1 

15.11 

13 . 44 

J 

15 . 16 

12.31 

1 

13.85 

8.94 

1 

11.23 


ORIGINAL PAGE rs 
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» rtarrr 


File: LINBEAR 2 . FRN 2,635 .a.. 


6 - 22-94 


Page 1 


49121.9 , 

49533.3 , 

49978.6 , 

50536.8 , 

51267.6 , 

5 1 886 .9 , 

52012.0 , 

52324.3 , 

52855.0 , 

53242.1 , 

54257.9 , 

54648.6 , 

54765.9 , 

55085 .9 , 

55427.0 , 

55797.6 , 

56398.3 , 
56782 . 1 . , 

56931.3 , 

57141.6 , 

57460.7 , 

57748.7 , 

58041.7 , 

58216.6 , 

58338.6 , 

58734.0 , 

58873.3 , 

59283.9 , 

59657 . 1 , 

60050 . 6 , 

60181.3 , 


0.0250 , 
0.0239 , 
0.0233 , 
0.0462 , 
0.0345 , 
0.0244 , 
0.0260 , 
0.0228 , 
0 . 0225 , 
0.0460 , 
0.0299 , 
0.0230 , 
0 . 0265 , 
0.0229 , 
0.0227 , 
0.0459 , 
0.0316 , 
0.0239 , 
0.0282 , 
0.0234 , 
0.0224 , 
0.0461 , 
O . 03 1 6 , 
0.0217 , 
0.0252 , 
0.0238 , 
0.0249 , 
0.0458 , 
0.0385 , 
0.0289 , 
0.0294 , 


0 . 0272 
0.0252 
0.0260 
0.0107 
0.01 89 
0.0247 
0.0268 
0.0271 
0.0275 
0.0114 
0 .0209 
0 . 0254 
0.0260 
O . 0270 
0.0274 
■ 0.0102 
0.0197 
0. 0245 
0.0280 
0 . 0265 
0.0273 
0.0108 
0 . 0203 
0 . 0257 
0.0281 
0 • 0267 
0. 0270 
0.0118 
0.0177 
0 . 0240 
0.0253 


0.0305 


0.0275 


0 


0 

9 

0 

9 

o 

600 

0 . 0280 


0. 0252 


61 


60 

9 

59 

9 

0.0271 


0.0228 


61 


61 

9 

59 

9 

1 900 

0. 0000 


0 . 0000 


601 


113 

9 

114 

9 

601 

0.01 55 


0.0438 


302 


298 

f 

299 

9 

600 

0.0207 


0.0289 


301 


298 

9 

298 

9 

1800 

0.0275 


0.0292 


0 


0 

9 

0 

9 

o 

0.0257 


0.0254 


61 


59 

9 

59 

9 

601 

0 . 0254 


0.0229 


60 


59 

9 

60 

9 

ISO ! 

0 . 0000 


0 . 0000 


602 


115 

9 

113 

1 

601 

0 .0166 


0 . 040 1 


303 


297 

9 

299 

9 

60 : 

0. 0205 


0.0278 


299 


299 

9 

299 

9 

1 79 > 

0.0276 


0.0301 


0 


0 

9 

0 

9 

i * 

0.0260 


0 . 0257 


61 


60 

9 

59 

9 

594 

0.0257 


0.0232 


60 


60 

9 

59 

9 

1 30 1 

0 . 0000 


0 . 0000 


599 


1 15 

9 

114 

9 

60 1 

0. 0164 


0.0416 


302 


300 

9 

299 

9 

601 

0.0206 


0 . 0288 


301 


301 

9 

300 

9 

1801 

0.0233 


0.0335 


0 


0 

9 

0 

9 

c 

5*9 

0. 0259 


0 . 0264 


61 


60 

9 

59 

9 

0. 0256 


0 . 0230 


61 


59 

9 

61 

9 

180 '.: 

0. 0000 


0 . 0000 


600 


1 15 

y 

1 15 

9 

602 

0.0159 


0.0415 


302 


300 

9 

300 

9 

600 

0 . 0208 


0.0262 


301 


300 

9 

298 

9 

1802 

0 . 0265 


0.0291 


0 


0 

9 

0 

9 

0 

0.0257 


0. 0267 


61 


60 

9 

59 

9 

60 1 

0.0254 


0.0256 


60 


60 

9 

59 

9 

1802 

0 . 0000 


0 . 0000 


602 


113 

9 

114 

9 

607 

0.0127 


0. 0486 


302 


300 

9 

299 

9 

5 * c 

0.0183 


0.0353 


302 


300 

9 

301 

9 

1307 

0.0243 


0 . 0350 


0 


0 

9 

0 

9 

f 
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LI 


(^ t oool H PEAP &AN E-* 

BEARING TEST NO 2 @ A ZONE 6-20-94 
(5 runs were performed on each case ) 
B1....B4 converted to INCH 
LI ...L4 in LBf 


TIME 

49121.9 

49533.3 
49978 6 
50536.8 

51267.6 
51886 9 
52012 0 

52324.3 

52855.0 

53242.1 
54257 9 

54648.6 
54765 9 
550859 
55427.0 

55797.6 
56398 3 
56782 1 
56931 3 
57141 6 

57460.7 
57748 7 

58041.7 
582166 
58338 6 
58734.0 
58973 3 
59283 9 
59657 1 
60050.6 

60181.3 


B1 

0 0250 
0.0239 
0.0233 
0 0462 
00345 
0.0244 
0 0260 
0 0228 
0.0225 
0.0460 
0.0299 
0.0230 
0 0265 
0.0229 
0.0227 
0.0459 
0 0316 
0 0239 
0 0282 
0.0234 
0 0224 
0.0461 
0.0316 
0 0217 
0 0252 
0.0238 
0 0249 
0.0458 
0 0385 
0 0289 
0.0294 


B2 

0 0272 
00252 
0.0260 
0.0107 
0.0189 
.0.0247 
0 0268 
0.0271 
00275 
0.0114 
0.0209 
0.0254 
0 0260 
0.0270 
00274 
0.0102 
0.0197 
0.0245 
0.0280 
0 0265 
0 0273 
0.0108 
0.0203 
0.0257 
00231 
0.0267 
0 0270 
0.0118 
0 0177 
0.0240 
0.0253 


B3 

B4 

0 0305 

0 0275 

0.0280 

0 0252 

0.0271 

0 0228 

0.0000 

0 0000 

0.0155 

0.0438 

0.0207 

0.0289 

00275 

0.0292 

0.0257 

0 0254 

0.0254 

0 0229 

00000 

0.0000 

0.0166 

0 0401 

00205 

0 0278 

00276 

0.0301 

0.0260 

0 0257 

0.0257 

0.0232 

0 0000 

0.0000 

0 0164 

0.0416 

00206 

0 0283 

00233 

0 0335 

0.0259 

0 0264 

00256 

0 0230 

0.0000 

0.0000 

0.0159 

0.0415 

0.0208 

0 0262 

0.0265 

0.0291 

0 0257 

0 0267 

0.0254 

0.0256 

0.0000 

0.0000 

0.0127 

0.0486 

0.0183 

0.0353 

0.0243 

0.0350 


10 

Do wti 


I — ^ c*. /Is 


LI 

L2 

L3 

L4 

CAGE 

0 

0 

0 

0 

0 

61 

60 

59 

600 

5 

61 

61 

59 

1800 

4 

601 

113 

114 

601 

6 

302 

298 

299 

600 

2 

301 

298 

298 

1800 

1 

0 

0 

0 

0 

0 

61 

59 

59 

601 

5 

60 

59 

60 

1802 

4 

602 

115 

113 

601 

6 

303 

297 

299 

600 

A 

< 

299 

299 

299 

1799 

1 

0 

0 

0 

0 

C 

61 

60 

59 

599 

5 

60 

60 

59 

1801 

4 

599 

115 

114 

601 

6 

302 

300 

299 

601 

2 

301 

301 

300 

1801 

1 

0 

0 

0 

0 

0 

61 

60 

59 

599 

5 

61 

59 

61 

1800 

4 

600 

115 

115 

602 

6 

302 

300 

300 

600 

2 

301 

300 

298 

1802 

1 

0 

0 

0 

0 

0 

61 

60 

59 

601 

5 

60 

60 

59 

1802 

4 

602 

113 

114 

603 

6 

302 

300 

299 

599 

2 

302 

300 

301 

1803 

1 

0 

0 

0 

0 

0 
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BENTLY PROBE CALIBRATION 

6 - 17-94 



(SllOA) 39V110A indino 

10-13 




BENTLY PROBE CALIBRATION 

6 - 20-94 



(snoA) 3DVHOA indino 
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ATTACHMENT 4a 

GAUSS CODE TO CALCULATE FEREBA INPUT LOADS 


h 


file BFORCE . PRC 


DHM 


7-16-94 


This GAUSS code transforms the forces applied to the bearing 
tester to FEREBA forces and moments at the inner 
for post -test correlation as part of the ROLLING ELEMENT 
BEARING MECHANICS contract with MSFC. 


INPUT 

*/ 


R = ZEROS (4,4) ; 

R [1 # 1] 

3= 

1.0; 

R [1, 2] 

« 

1.0; 

R[l» 3] 

a 

1.0; 

R [2 , 4] 


1.0; 

R [3 , 1] 

3C 

3.600; 

R [3 , 2] 

« 

-1.800; 

R [3 , 3] 

* 

-1.800; 

R [3 , 4] 


-0.15; 

R [4 , 2] 

- 

-3.118; 

R [4 , 3 ] 


3.118; 

L = ZEROS (4,1) ; 

L [1 , 1] 

- 

608.8; /* 

L [2 , 1] 

* 

122.0; /* 

L [3 , 1] 

= 

122.2; /* 

L [4 , 1 ] 

a 

601.6; /* 


AVG LI FORCE FOR 
AVG L3 FORCE FOR 
AVG L2 FORCE FOR 
AVG L4 FORCE FOR 


CASE 

6 

TEST 

2 

*/ 

CASE 

6 

TEST 

2 

*/ 

CASE 

6 

TEST 

2 

*/ 

CASE 

6 

TEST 

2 

*/ 


/* 

CALCULATE FEREBA FORCES AND MOMENTS 

*/ 

FFM * R*L; 

/* 

OUTPUT DATA 

*/ 

OUTPUT FILE = BFORCE. OUT; 
OUTPUT RESET; 

FORMAT /Ml /RD 12,4; 


PRINT FFM; 
PRINT; 
PRINT; 
PRINT; 

PRINT R; 
PRINT; 

PRINT L; 
OUTPUT OFF; 
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ATTACHMENT 4b 

CALCULATED FEREBA INPUT LOADS FROM GAUSS CODE 



FILE BFORCE . OUT 

DHM 

7-18-94 

— 

FEREBA APPLIED 

FORCES AND 

MOMENTS FOR TEST 2 (CASE 

- 

FFM 




853.0000' 
601 . 6000 




1661.8800 




0.6236 



- 

R 

1.0000 ■ 
0.0000 
3.6000 
0.0000 

1.0000 

0.0000 

-1.8000 

-3.1180 

1.0000 0.0000 

0.0000 1.0000 

-1.8000 -0.1500 

3.1180 0.0000 


L 

608.8000 

122.0000 

122.2000 

601.6000 
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ATTACHMENT 5a 


GAUSS CODE TO CALCULATE "MEASURED" SHAPT DISPLACEMENTS 


file BDEFLA1 . PRC 


DHM 


7-16-94 


This GAUSS code calculates the inner ring/shaft deflections 
in FEREBA coordinates from positions measured in the bearing 
static test performed as part of the ROLLING ELEMENT BEARING 
MECHANICS contract with MSFC. 

This program is hardwired for load case 1 of test A without 
deadband. 


INPUT 

*/ 


ZEROS (4,4) ; /* MRDV - R*FDEFL */ 


R[l,l] 

* 

1.0; 

R [2 , 1] 

* 

1.0; 

R [3 , 1] 

* 

1.0; 

R [4 , 2] 

* 

1.0; 

R [1 , 3] 

« 

- 1.993; 

R [2 , 3 ] 

* 

0.997; 

R [3 , 3] 

* 

0.997; 

R [4 , 3] 

* 

- 2.10; 

R [2 , 4] 

* 

- 1.726; 

R [3 , 4] 

« 

1.726; 

RINV 


INV(R) ; 


B = ZEROS (4,5) ; 
LOAD CASE 1 


/* LOADED POSITIONS FOR 5 TRIALS IN 
IN FEREBA COORDINATES *f 

'/ 


'/ 


B [ 1 , 1] 

sc 

0.0197 

/* 

POSITION AT GAGE 

B1 

B [1, 2] 

= 

0.0198 




B [1,3] 


0.0197 




B [ 1 , 4 ] 

■ 

0.0196 




B [1 , 5] 

m 

0.0201 




B [2 , 1] 

« 

0.0317 

/* 

POSITION AT GAGE 

B2 

B [ 2 , 2 ] 


0.0316 




B [2 , 3] 


0.0316 

/* 

GAGE B3 IN TEST 

*/ 

B [2 , 4 ] 

> 

0.0317 




B [2 , 5] 

S 

0.0312 




B [3 , 1] 

s 

0.0272, 

/* 

POSITION AT GAGE 

B3 

B [3 , 2] 

B 

0.0272 




B [3 , 3 ] 

* 

0.0273 

/* 

GAGE B2 IN TEST 

*/ 

B [3 , 4] 

* 

0.0273 




B [3 , 5] 


0.0271 





'/ 
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B [4 , 1] 
B [4, 2] 
B [4 , 3 ] 
B [4 , 4 ] 
B [4 , 5] 


0.0238 

0.0240 

0.0239 

0.0238 

0.0244 


/* POSITION AT GAGB B4 */ 


BO = ZEROS (4, 5); /* UNLOADED POSITIONS FOR 


B0 [1,1] 

as 

0.0292; 


/* 

POSITION 

AT 

GAGE 

B1 

B0 [1,2] 


0.0291; 







BO [1,3] 

«= 

0.02918; 






B0 [1,4] 

= 

0.0292 







B0 [1,5] 

= 

0.0292 







B0 [2,1] 

s 

0.0313 


/* 

POSITION 

AT 

GAGE 

B2 

B0 [2 , 2] 

tm 

0.0317 







B0 [2,3] 

= 

0.03155; 

/* 

GAGE B3 

EN TEST 

*/ 

B0 [2,4] 

m 

0.0316 







B0 [2,5] 

s 

0.0316 







B0 [3 , 1] 

SB 

0.0280 


/* 

POSITION 

AT 

GAGE 

B3 

BO [3,2] 

* 

0.0281 







BO [3,3] 


0.02805; 

/* 

GAGE B2 

IN TEST 

*/ 

BO [3,4] 

« 

0.0281 







BO [3,5] 

JZ 

0.0280 







BO [4,1] 

JE 

0.0318 


/* 

POSITION 

AT 

GAGE 

B4 

BO [4,2] 

« 

0.0315 







BO [4,3] 

SB 

0.03165; 






BO [4,4] 

= 

0.0316; 






BO [4,5] 

= 

0.0317 

/ 







5 TRIALS 

*/ 

*/ 




'/ 


CALCULATE MEAN RELATIVE DEFLECTION VECTOR FOR THIS LOAD CASE 

*/ 

MRDV = ZEROS (4,1); 

i - l; 

DO WHILE i LE 4; 

j - i; 

SUM = 0.0; 

DO WHILE j LE 5; 

SUM = SUM + B0 [i , j ] - B [ i , j ] ; 

j - j + i; 

ENDO; 

MRDV [ i , 1 ] - SUM*0 . 2 ; 
i » i + 1; 

ENDO; 

/* 

CALCULATE FEREBA DEFLECTIONS FOR THIS LOAD CASE 

*/ 

FDEFL - RINV*MRDV; 
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OUTPUT DATA FOR THIS LOAD CASE 

*/ 

OUTPUT FILE = BDEFLAl.OUT; 
OUTPUT RESET; 

FORMAT /Ml /RD 12,4; 

PRINT FDEFL; 

PRINT; 

PRINT; 

PRINT B; 

PRINT; 

PRINT MRDV ; 

PRINT; 

PRINT R; 

PRINT; 

PRINT BO ; 

PRINT; 

OUTPUT OFF; 
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ATTACHMENT 5b 

CALCULATED "MEASURED" SHAFT DISPLACEMENTS FROM GAUSS CODE 


FILE BDEFLA1 . OUT DHM 7-18-94 

"MEASURED" SHAFT DEFLECTIONS FOR TEST 1 (CASE 1) 


FDEFL 

0.0034 

0.0014 

-0.0030 

0.0002 





B 

0.0197 

0.0317 

0.0272 

0.0238 

0.0198 

0.0316 

0.0272 

0.0240 

0.0197 

0.0316 

0.0273 

0.0239 

0.0196 

0.0317 

0.0273 

0.0238 

0.0201 

0.0312 

0.0271 

0.0244 

MRDV 

0.0094 

- 0.0000 

0.0008 

0.0077 





R 

1.0000 

1.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

-1.9930 

0.9970 

0.9970 

-2.1000 

0.0000 

-1.7260 

1.7260 

0.0000 


B0 

0.0292 

0.0313 

0.0280 

0.0318 

0.0291 

0.0317 

0.0281 

0.0315 

0.0292 

0.0316 

0.0280 

0.0316 

0.0292 

0.0316 

0.0281 

0.0316 

0.0292 

0.0316 

0.0280 

0.0317 


TABLE 1 


TRANSDUCER LOCATIONS FOR BEARING TEST RIG 


TRANSDUCER 

DIRECTION 

XF LOCATION 

YF LOCATION 

ZF LOCATION 

LI 

Axial 

N/A 

0.000 

3.600 

L2 

Axial 

N/A 

3.118 

-1.800 

L3 

Axial 

N/A 

-3.118 

-1.800 

L4 

Radial 

0.15 

N/A 

N/A 

Bl 

Axial 

N/A 

0.000 

-1.993 

B2 

Axial 

N/A 

1.726 

0.997 

B3 

Axial 

N/A 

-1.726 

0.997 

B4 

Radial 

2.10 

N/A 

N/A 


Notes : 

(1) Li refers to load cell while Bi refers to deflection transducer. 

(2) The origin of the FEREBA coordinate system is at bearing center. 

(3) YF and ZF coordinates obtained from drawings and checked on rig. 

(4) XF coordinates measured from rig and subject to +/- 0.10" error. 
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TABLE 2 

TEST RIG AND FEREBA APPLIED LOADS 

TEST NO. 1 (NO DEADBAND) 
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(2) Fx is axial; Fy and Fz are radial; My is moment about FEREBA Y axis; Mz is about Z. 

(3) The requested radial force for Case. 3 (3000 lb) could not be achieved by the test rig 



TEST RIG AND FEREBA APPLIED LOADS 

TEST NO. 2 (O.OOOI ” DEADBAND) 
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TABLE 4 

MEASURED" AND FEREBA SHAFT DISPLACEMENTS 

TEST NO. 1 (NO DEADBAND) 
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(2) Shaft displacements, calculated by FEREBA from equivalent forces applied to test ri 
less the axial shaft displacement of 0.0006 inch due to the 24 lb preload. Aster is 
indicate that FEREBA did not converge for these load cases. 





TABLE 5 

MEASURED" AND FEREBA SHAFT DISPLACEMENTS 

TEST NO. 2 {0.0001" DEADBAND) 
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